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Abstract 

In  this  research  the  close  formation  flight  control  problem  is  addressed.  The  forma¬ 
tion  consists  of  a  lead  and  wing  aircraft,  where  the  wing  flies  in  close  formation  with  the 
lead,  such  that  the  lead’s  vortices  produce  aerodynamic  coupling  effects,  and  a  reduction  in 
the  formation’s  drag  is  achieved.  A  controller,  i.e.,  a  formation- hold  autopilot  for  the  wing 
aircraft,  is  designed  such  that  the  formation’s  geometry  is  maintained  in  the  face  of  lead 
aircraft  maneuvers.  In  the  formation  flight  control  system,  the  wing  and  lead  aircraft  dy¬ 
namics  are  coupled  due  to  kinematic  effects,  and,  in  the  case  of  close  formations,  additional 
aerodynamic  coupling  effects  are  introduced.  In  the  research  these  additional  aerodynamic 
coupling  effects  are  properly  modeled.  The  most  significant  aerodynamic  coupling  effect 
introduced  by  close  formation  flight  entails  the  coupling  of  the  lateral/ directional  channel 
into  the  altitude-hold  autopilot  channel.  It  is  shown  that  formation  hold  autopilots  de¬ 
signed  ignoring  the  aerodynamic  coupling  effects,  yield  satisfactory  performance  in  close 
formation  flight. 
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CLOSE  FORMATION  FLIGHT  CONTROL 

/.  Introduction 

1.1  Overview  of  Thesis 

The  introduction  and  review  of  previous  research  related  to  this  subject  are  discussed 
in  Chapter  L  The  equations  of  motion  describing  the  kinematics  between  the  aircraft  in 
the  formation  as  well  as  the  flight  control  system  equations  for  the  aircraft  are  provided  in 
Chapter  IL  The  derivation  of  the  Lead  aircraft’s  effects  on  the  Wing  aircraft  is  developed 
in  Chapter  III.  This  includes  derivations  for  upwash  and  sidewash  effects  caused  by  the 
Lead’s  wing  vortex  on  the  Wing  aircraft  and  the  development  of  the  modifled  Wing  aircraft 
flight  control  system.  In  Chapter  IV  the  close  formation  flight  controller  along  with  the 
complete  simulation  model  is  developed.  An  analysis  of  stability  and  robustness  is  pre¬ 
sented.  Simulations  are  conducted  and  a  performance  evaluation  of  results  is  presented  in 
Chapter  V.  Analysis,  conclusions,  and  recommendations  for  further  research  are  presented 
in  Chapter  VI. 

1.2  Background 

Previous  research  has  been  performed  by  several  AFIT  graduate  students  on  the 
subject  of  formation  flight  control  [1]  -  [12].  The  motivation  for  their  work  was  based 
primarily  on  the  needs  of  special  operations  forces  (SOF).  SOF  missions  require  pilots  to 
fly  many  hours  over  long  distances  in  a  large  formation  usually  in  the  dark  at  low  altitudes 
and  at  times  in  poor  weather.  These  flying  requirements  combined  with  the  complexities 
of  todays  modern  aircraft  are  very  stressful  and  can  lead  to  pilot  saturation.  A  second 
motivation  for  this  research  was  the  need  to  fly  multiple  unmanned  aerial  vehicles  (UAV) 
in  a  formation.  These  UAVs  may  be  flying  a  surveillance  mission  or  patroling  a  no-fly  zone. 
Current  UAVs  are  limited  in  the  number  of  payloads  they  can  carry  at  one  time.  This 
makes  it  necessary  to  fly  multiple  UAVs  to  obtain  a  complete  set  of  data.  A  formation 
controller  offers  two  potential  benefits.  First,  if  used  for  close  formations,  it  can  reduce  fuel 
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consumption  of  the  formation;  secondly,  it  can  keep  the  aircraft  in  a  known  position  with 
respect  to  each  other,  making  correlation  of  data  easier  for  fusion  from  multiple  sensor 
payloads  in  multiple  UAV s.  In  general,  this  could  also  help  any  formation  of  aircraft  flying 
long  distances. 

In  previous  work,  the  formation  flight  control  problem  was  introduced  and  attention 
was  given  to  the  kinematic  coupling  effects-see,  e.g.,  [1]-[12].  Multiple  versions  of  formation 
flight  controllers  were  designed  and  developed.  These  controllers  consisted  of  both  feedback 
and  feedforward  designs  as  well  as  energy  conserving  designs.  Simulations  demonstrated 
that  these  controllers  work  well  for  large  formations  in  which  the  separations  are  more  than 
a  wingspan.  However,  they  were  not  designed  for  close  formations,  in  which  the  lateral 
separation  between  aircraft  is  less  than  a  wingspan. 

1,3  Problem  Statement 

An  aircraft  while  flying  causes  spiral  like  disturbances  in  the  air.  Two  disturbances 
are  projected  backward  perpendicular  to  it’s  wing  and  one  is  parallel  to  the  wing.  The 
spiral  disturbances  are  cylindrical  in  shape  and  dissipate  slowly  over  time.  The  air  in 
the  disturbances  flows  in  a  circular  pattern  starting  at  the  tip  of  the  wing,  flowing  up 
circularly  towards  the  body  of  the  aircraft  and  continuing  to  circularly  flow  down  back  to 
the  wing  tip.  This  leads  to  a  clockwise  flow  in  the  disturbance  created  by  the  left  wing 
and  a  counter-clockwise  flow  from  the  disturbance  created  by  the  right  wing,  as  observed 
from  behind.  These  disturbances  created  by  the  wing  of  the  aircraft  are  often  called  a 
horseshoe  vortex,  see  Figure  1.1.  The  air  flow  directly  between  the  two  disturbances  is 
almost  purely  downwash.  The  flow  of  the  air  on  the  outside  of  the  horseshoe  vortex  is  a 
mixture  of  up  wash  and  si  dewash.  Aircraft  flying  into  a  vortex  created  by  another  aircraft 
can  experience  substantial  aerodynamic  forces.  The  FAA  has  even  developed  a  separation 
criteria  to  reduce  the  danger  of  accidents  from  other  aircraft  flying  into  the  wake  caused 
by  aircraft  landing  and  taking  off  [16]. 

Prandtl,  a  German  scientist  was  one  of  the  first  to  analyze  the  effects  of  these  vortices 
on  other  aircraft.  It  has  been  hypothesized  that  geese  fly  in  their  inverted  ‘V’  formation 
to  take  advantage  of  the  upwash  created  by  the  vortex  from  each  leading  goose.  Another 
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aircraft  flying  in  the  vortex  can,  if  positioned  correctly,  experience  a  reduction  in  drag  and 
an  increase  in  lift  [13].  Therefore,  alternating  aircraft  in  and  out  of  the  lead  position  can 
potentially  increase  the  range  and  endurance  of  a  formation  of  aircraft.  The  formation  of 
the  two  aircraft  is  shown  in  Figure  1.2.  The  aircraft  creating  the  vortex  is  the  Lead  and 
the  aircraft  flying  behind  and  into  the  vortex  is  the  Wing.  Indeed,  the  vortex  shed  by 
the  Lead  aircraft  induces  an  upwash  on  the  Wing  aircraft’s  wing,  which  is  responsible  for 
reducing  the  Wing’s  induced  drag. 

The  difficulty  in  taking  advantage  of  this  upwash  for  drag  reduction  by  another 
aircraft  is  the  fact  it  requires  very  precise  positioning  with  respect  to  the  aircraft  generating 
the  vortex.  There  is  a  small  cylindrical  window  with  a  radius  of  about  ten  percent  of  the 
wingspan  in  which  a  benefit  can  be  achieved  for  the  Wing  aircraft  [13].  However,  these 
vortices  have  substantial  force  and  the  flight  control  system  needs  to  be  able  to  reject  these 
disturbances  and  keep  the  aircraft  stable  and  in  the  proper  location. 


b 


Figure  1.1  Horse  Shoe  Vortex 


This  thesis  is  a  change  in  direction  from  the  previous  work  on  formation  flight  control. 
The  main  problem  focus  in  this  research  is  to  analyze  the  effects  of  close  formation  flying 
on  the  wing  aircraft’s  flight  control  system  and  develop  a  close  formation  flight  controller 
robust  enough  to  handle  these  effects  and  maintain  the  close  formation.  This  will  primarily 
include  the  need  to  analyze  the  vortex  created  by  the  Lead  and  determine  its  change  in 
strength  based  on  a  change  in  position  in  the  formation.  Secondly,  the  necessary  modi¬ 
fications  to  be  applied  to  a  formation  flight  control  system  to  maintain  a  close  formation 
during  flight  must  be  determined. 

1.4  Research  Objectives  and  Questions 

The  primary  objective  of  this  research  is  to  apply  a  formation  flight  controller,  previ¬ 
ously  developed  for  large  formation,  to  this  problem  in  which  a  close  formation  is  required. 
Flight  control  system  models  previously  developed  will  also  be  modified  for  use  in  this 
research. 

The  primary  question  to  be  answered  from  this  work  is,  “can  a  close  formation  flight 
controller  developed  without  due  consideration  of  vortex  effects  control  a  close  forma¬ 
tion?”.  A  secondary  question  is:  “which  is  stronger,  kinematic  cross  coupling  or  vortex 
aerodynamic  cross  coupling?”. 

1.5  Assumptions 

There  are  ten  main  assumptions  listed  in  order  of  development.  These  assumptions 
reduce  the  complexity  and  scope  of  the  problem.  This  allows  a  focused  study  on  the 
research  objectives.  These  assumptions  do  not  negate  the  reasonableness  of  the  results 
achieved.  An  unlisted  objective  is  to  keep  this  research  as  realistic  as  possible.  Hopefully, 
the  results  can  be  applied  to  the  development  of  an  actual  formation  flight  controller  for 
use  in  future  aircraft  and  UAVs. 

•  Formation  flight  controller  resides  on  Wing  aircraft. 

•  Each  aircraft  is  modeled  as  a  point  mass.  This  means  that  moment  effects  are  not 

considered. 
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WfNG 

Figure  1.2  Formation  Geometry 

•  The  aircraft  flight  control  system  is  assumed  to  be  composed  of  three  independent 
autopilots:  Mach  hold,  heading  hold,  and  altitude  hold. 

•  Rotating  coordinate  frame  is  attached  to  the  wing  aircraft:  x  out  the  nose,  y  out  the 
right  wing,  and  z  down,  see  Figure  1.2. 

•  A  sensor  capable  of  measuring  the  separation  distances  as  well  as  the  velocity,  heading 
and  altitude  of  the  lead  aircraft  is  attached  to  the  wing  aircraft. 

•  The  Horseshoe  vortex  is  modeled  as  two  filaments  originating  near  each  wing  tip  and 
proceeding  backward  for  an  infinite  distance.  The  magnitude  of  the  parallel  filament 
effect,  that  normally  completes  the  horseshoe,  is  assumed  to  be  negligible.  These 
filaments  are  shown  in  Figure  1.1. 

•  The  disturbances  created  by  the  Wing  on  the  Lead  are  ignored. 

•  The  Wing  aircraft  is  assumed  to  be  in  the  most  “optimal”  position.  “Optimal”  means 
that  the  Wing  is  in  a  position  of  maximum  reduction  in  drag. 
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•  The  wing  flight  control  system  is  trimmed  for  straight  and  level  flight  at  this“optimal” 
position  in  the  formation. 

•  The  Lead  maneuvering  envelope  is  :  +/-  30  degree  heading  change,  +/-  50  ft/sec 
velocity  change,  and  +/-  1000  ft  altitude  change. 

1.6  Approach/ Methodology 

Systems  developed  by  previous  AFIT  students  axe  used  as  the  building  blocks  for  this 
research  [11],  [1],  and  [3].  These  blocks  include  a  Flight  Control  System  (FCS)  and  Forma¬ 
tion  Flight  Controller  (FFC).  The  FCS  model  originally  developed  by  Rohs  [11]  was  a  first 
order  model  based  on  C-130  aircraft  models  provided  by  Lockheed.  Buzogany  [1]  developed 
a  second  order  flight  control  system  model  through  system  identification  techniques.  This 
second  order  model  was  reduced  from  a  fourteenth  order  model  that  was  developed  based 
on  actual  aircraft  specifications  provided  by  Lockheed.  The  complete  flight  control  system 
used  in  this  research  is  composed  of  the  first  order  Heading  hold  autopilot  developed  by 
Rohs  and  the  first  order  Mach  hold  autopilot  and  the  second  order  Altitude  hold  model  de¬ 
veloped  by  Buzogany.  These  three  autopilots  constitute  the  complete  flight  control  system 
used  for  both  the  Lead  and  Wing  aircraft.  In  his  thesis  [3],  Dargan  developed  a  FFC  using 
mixed  feedback  compensation.  The  controller  was  a  three  channel  controller  for  controlling 
the  horizontal  x  and  y  separation  distances  and  the  vertical  z  separation  distance  in  the 
formation.  The  formation  controller  contained  a  proportional  plus  integral  (PI)  compen¬ 
sator  on  each  channel  with  a  linear  mixer  for  the  inputs  to  the  PI  compensators.  The 
flight  control  system  developed  by  Rohs  and  Buzogany  and  the  formation  flight  controller 
developed  by  Dargan  are  modified  and  used  as  the  foundations  for  this  research. 

The  research  analysis  begins  with  the  development  of  a  simulation  using  the  first  order 
flight  control  system  developed  by  Rohs  and  the  formation  flight  controller  developed  by 
Dargan.  An  initial  simulation  model  is  developed  in  Simulink  [17].  This  model  is  used 
to  verify  the  results  achieved  by  Dargan.  This  ensures  the  flight  control  system  model 
and  the  formation  flight  controller  are  accurately  modeled.  Once  confidence  is  gained  in 
the  complete  simulation  model,  the  FCS  models  are  adjusted  to  represent  an  F-16  class 
aircraft.  The  derivation  for  the  effects  of  the  horseshoe  vortex  are  developed  for  a  close 
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formation.  These  effects  are  then  integrated  into  the  simulation  as  disturbances  to  the 
wing  and  the  FFC  is  assessed.  The  FFC  is  adjusted  to  obtain  the  best  performance. 

Simulations  are  performed  using  Matlab  version  4.21  Simulink  package  by  The  Math- 
Works  [17].  Simulink  allows  use  of  simple  blocks  to  form  a  frequency  domain  representation 
of  a  complex  system.  The  non-linear  differential  equations  are  integrated  using  a  fourth 
order  Runge-Kutta  technique  and  an  Adams  technique  contained  in  Matlab  Simulink. 
Both  have  an  error  limiting  capability  and  step  size  control  to  ensure  accurate  results  of 
integrations. 
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IL  Model  and  Simulation  Development 


This  chapter  outlines  the  models,  kinematic  equations,  and  basic  controller  used  in 
the  nonlinear  simulation.  The  majority  of  Section  2.1  has  been  taken  from  [12].  Models 
and  figures  have  been  changed  to  refiect  the  use  of  an  F-16  class  aircraft  instead  of  a  C-130 
class  aircraft.  Sections  2.2,  2.3  and  2.4  have  been  taken  directly  from  [12];  however,  some 
figures  have  been  modified  slightly  for  the  reasons  stated  above  and  a  few  paragraphs  have 
been  added  to  provide  additional  information.  Section  2.5  contains  a  description  of  the 
controller  developed  by  Dargan  and  walks  through  his  controller  design  process.  Sufficient 
information  is  provided  to  enable  the  reader  to  reproduce  the  results  of  this  research. 

2J  Aircraft/ Autopilot  Models 

The  development  of  the  aircraft/autopilot  models  was  accomplished  by  Rohs  [11]  and 
Buzogany  [1].  Since  autopilot  models  were  not  explicitly  available,  models  were  obtained  by 
designing  a  custom  autopilot  system  around  C-130H  aircraft  models  provided  by  Lockheed. 
The  autopilot  system  included  Heading-hold,  Mach-hold,  and  Altitude-hold  autopilots. 
Due  to  the  inherent  properties  of  aircraft  autopilots,  the  models  exhibited  overdamped, 
decoupled  responses  [11].  Using  system  identification  techniques,  first-order  models  were 
developed  using  rate  limited,  decoupled,  first-order  differential  equations.  The  first-order 
aircraft /autopilot  models  are  specified  as: 


(2.1) 
(2.2) 
(2.3) 

where, 

=  heading  time  constant 
Tt,  =  velocity  time  constant 
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Th  =  altitude  time  constant 

The  first-order  aircraft/autopilot  models  are  illustrated  in  Figure  2.1. 


Figure  2.1  First-Order  Aircraft  Models 


Second-order  aircraft/autopilot  models  were  developed  by  Buzogany  [1]  in  order  to 
more  accurately  represent  the  “true”  aircraft /autopilot  system.  The  most  disturbing  prob¬ 
lem  with  first-order  models  is  the  presence  of  instantaneous  heading  rate  and  vertical  veloc¬ 
ity  changes.  Figure  2.2  compares  an  overdamped  second-order  response  with  a  first-order 
response.  While  heading  and  altitude  response  was  significantly  improved  using  second- 
order  models,  it  was  found  that  velocity  response  was  modeled  more  precisely  using  a 
first-order  model  with  a  larger  time  constant.  Thus,  the  “second-order”  models  and  time 
constants  are  specified  as: 


V  = 

^  = 
h  = 


-Iv  +  iy. 

T~v  '^v 

-  f  —  +  — )  i’ - - — Ip  +  — - — V’c 

-  -t-  — )  h  -  -^h  -t- 


(2.4) 

(2.5) 

(2.6) 
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Figure  2.2  Comparison  of  First  and  Second-Order  Responses 


There  are  multiple  ways  to  model  a  second  order  system  containing  nonlinear  sat¬ 
urations  caused  by  actuators.  The  models  illustrated  in  Figure  2.3  were  found  to  be  the 
best  representation  of  a  true  aircraft  response.  Comparison  of  these  models  with  models 
not  including  saturations  are  almost  identical  for  the  envelope  of  maneuvers  analyzed  in 
this  research. 

The  complete  flight  control  system  is  composed  of  Rohs’  first  order  Heading  hold 
autopilot  model  and  Buzogany’s  first  order  Mach-hold  and  second  order  Altitude-hold 
autopilot  models. 


V  =  --V+-V,  (2.7) 

'^V 

i’  =  V’+— V’c  (2.8) 

h  =  -(l-  +  ±)h-—h+-^K  (2.9) 

These  models  represented  a  C-130  class  aircraft.  The  C-130  aircraft  is  a  transport  and 
responds  slowly  in  comparison  to  an  F-16  class  aircraft.  Since  Rohs’  and  Buzogany’s 
models  are  used  in  this  research,  the  time  constants  and  saturation  values  need  to  reflect  the 
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Figure  2.3  Second-Order  Aircraft /Autopilot  Models 
Table  2.1  F-16  Class  Composite  Aircraft /Autopilot  Parameters 


Parameter 

Value 

0.33  sec 

Tv 

5  sec 

0,307  sec 

n. 

3.843  sec 

performance  characteristics  for  an  F-16  class  aircraft  versus  a  C-130.  First,  this  consisted 
of  simply  halving  the  time  constants  of  the  C-130  models  to  speed  up  the  system  response. 
This  should  be  reasonable  since  an  F-16  class  aircraft  responds  much  faster  than  C-130 
aircraft.  Secondly,  the  saturations  were  changed  to  match  the  maneuvering  capability  of 
an  F-16  class  aircraft.  In  general  the  saturations  are  about  two  to  three  times  greater  than 
those  of  a  C-130  class  aircraft.  These  new  saturation  values  should  also  be  reasonable  since 
the  maneuvering  envelope  of  an  F-16  is  much  greater  than  that  of  a  C-130.  Tables  2.1  and 
2.2  show  the  F-16  class  aircraft  time  constants  and  saturation  values  used  in  this  research 
respectively.  Figure  2.4  shows  the  complete  F-16  class  aircraft  flight  control  system  model 
used  in  this  research. 
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Figure  2.4  F-16  Class  Composite  Flight  Control  System  Model 


2,2  Equation  of  Coriolis 

In  order  to  express  the  formation  flight  control  problem  kinematics,  reference  frames 
must  be  established  in  the  inertial  and  wing  aircraft  frames.  Figure  2.5  shows  the  inertial 
and  wing  aircraft  frames  and  the  angular  rotation  vector.  The  subsequent  rotations  caused 
by  maneuvering  aircraft  within  the  formation  require  a  mathematical  function  relating 
vectors  in  different,  rotating  frames.  The  Equation  of  Coriolis  performs  this  function. 


Figure  2.5  Inertial  and  Rotating  Frames  of  Reference  [3:Figure  1.7] 

According  to  Blakelock  [15],  “The  motion  of  an  object  as  viewed  from  a  reference 
frame  is  equal  to  the  motion  as  seen  from  the  moving  frame,  plus  the  motion  resulting  from 
the  relative  angular  velocity  of  the  moving  frame  with  respect  to  the  reference  frame”  [15]. 
This  equation  provides  the  basis  for  converting  individual  aircraft  heading  and  velocity 
into  separation  distances.  This  is  an  essential  part  of  the  design  process.  The  Equation  of 
Coriolis  is: 

Ri  =  Rp+UipXRp  (2.10) 

where 
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•  Ri  =  the  vector  velocity  of  the  point  in  the  i  reference  frame 

•  JSp  =  the  vector  position  of  the  point  in  question  as  seen  from  the  p  reference  frame 

•  “  the  vector  angular  velocity  of  the  p  reference  frame  with  respect  to  the  i 
reference  frame 

•  Rp  =  the  vector  position  of  a  point  in  the  p  frame 
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2.3  Formation  Coordinate  System 

The  coordinate  system  used  in  this  thesis  is  identical  to  that  of  Dargan  [3],  Buzogany 
[1],  and  Reyna  [9].  The  analysis  of  the  system  kinematics  uses  two  coordinate  frames: 

•  Inertial  base  frame 

•  Rotating  reference  frame  centered  on  wing  aircraft 

The  base  frame  is  an  inertial  North-East-Down  system.  For  the  purposes  of  this 
research,  the  earth  is  considered  to  be  flat  and  inertially  flxed.  The  inertial  reference 
frame  and  separation  distances  are  shown  in  Figure  2.6. 


Figure  2.6  Inertial  Reference  Frame  and  Sepajation  Distances 

The  wing  aircraft  frame  is  centered  on  the  wing  aircraft.  The  x  axis  is  in  the  flight 
direction  (i.e.,  aligned  with  the  velocity  vector),  the  y  axis  points  out  the  starboard  wing, 
and  the  z  axis  points  toward  the  earth.  The  x  and  y  separation  distances  are  measured  in 
the  wing  frame  Figure  2.7. 
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2.4  Kinematic  Equations 


In  order  to  simulate  the  kinematics  associated  with  the  formation  flight  control  model, 
kinematic  equations  must  be  derived.  This  has  already  been  done  by  Dargan  [3] ,  Buzogany 
[1],  and  Reyna  [9].  Reyna’s  derivation  is  repeated  here  to  ensure  clarity  (from  [9],  pages 
3:6-10). 

Using  the  Equation  of  Coriolis  (Chapter  II),  the  velocity  of  the  lead  with  respect  to 
the  wing  has  been  found  by  Dargan  as 

V^L  =  V^ +u;^xK  (2.11) 

where  the  following  convention  is  followed: 

•  The  superscript  indicates  the  reference  frame. 

•  The  subscript  indicates  the  parameter  described  by  the  vector  or  a  relation  between 
two  parameters. 

For  example, 

•  ^WL  =  velocity  of  the  lead  aircraft  with  respect  to  the  Wing,  in  the  Wing’s  reference 
frame 

•  =  angular  velocity  of  the  wing  aircraft  in  the  wing’s  reference  frame 

•  ^WL  “  position  of  the  Lead  aircraft  with  respect  to  the  Wing  in  the  Wing’s  reference 
frame 

•  “  inertial  velocity  of  the  Wing  aircraft  in  its  own  reference  frame 

•  ==  inertial  velocity  of  the  Lead  aircraft  in  the  Wing’s  reference  frame 

The  development  of  the  kinematic  equations  is  based  on  the  geometry  defined  in 
Figure  2.8. 

The  following  relationships  are  defined: 
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where, 

'ipE  =  heading  error 

V[  =  velocity  of  the  lead  aircraft  in  its  own  reference  frame 

In  order  to  solve  Equation  (2.11),  F/  must  be  transformed  into  the  wing  reference 
frame.  A  Direction  Cosine  Matrix  (DCM),  developed  by  Dargan  [3],  is  used  to  perform 
the  rotation.  Dargan  found  the  DCM  to  be 


cos^E  -simj^E  0 
sin'ijjE  cos'ijjE  0 

0  0  1 


(2.18) 


is  found  using  by  substituting  (2.16)  and  (2.18)  into  the  following  equation 


=  cYv^  = 


ViCOSIpE 

VEsintpE 


0 


(2.19) 


Substituting  (2.13)-(2.15),  (2.17),  and  (2.19)  into  equation  (2.11)  yields 


- 

^WL  — 


ViCOStpE 

0 

‘  ' 

Vw 

0 

0 

VLsin'ij^E 

- 

0 

X 

yW 

- 

0 

+ 

0 

X 

0 

0 

_ 

0 

.  . 

0 

ViCOSIpE 

VisintpE 

- 

- 

0 

0 

(2.20) 


(2.21) 


Separating  (2.21)  into  scalar  components  yields 


=  VlCosiPb  +  -  Vw 


(2.22) 
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(2.23) 

(2.24) 


=  VisinipE  - 

=  0 

Equations  (2.22)  -  (2.24)  describe  the  kinematics  of  the  formation  (x  and  y  separations) 
in  terms  of  the  individual  aircraft’s  heading  and  velocity.  The  altitude  (z)  separations  are 
non-dynamic  and  are  simply  the  difference  between  the  aircraft  altitudes. 
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2,5  Formation  Flight  Controller 

Each  aircraft  is  controlled  through  the  inputs  to  the  flight  control  system.  These 
inputs  are  the  settings  normally  supplied  by  the  pilot.  The  control  strategy  implemented 
by  Dargan  was  a  two  tiered  structure  as  shown  in  Figure  2.9.  The  upper  tier  is  for  control 
of  the  formation  as  a  whole,  and  the  lower  tier  is  for  control  of  the  individual  aircraft 
within  the  formation  [3].  The  upper  tier  is  used  to  control  the  guidance  of  the  formation. 
This  consists  of  the  velocity,  heading,  and  altitude  for  the  formation.  The  lower  tier  is 
used  to  control  the  actual  formation  geometry.  This  consists  of  the  x,  y  and,  z  separation 
distances  between  the  Lead  and  Wing  aircraft  in  the  formation. 


Figure  2.9  FFC  Structure 


Dargan  [3]  separated  the  controller  into  three  control  channels:  longitudinal,  lateral 
and  vertical.  The  longitudinal  or  x-channel  was  composed  of  the  longitudinal  x  separation 
distance  and  the  relative  velocity  between  the  lead  and  wing  aircraft.  The  lateral  or  y- 
channel  was  composed  of  the  lateral  y  separation  distance  and  the  relative  heading  angle 
between  the  lead  and  wing  aircraft.  The  vertical  or  z-channel  was  simply  the  separation 
distance  between  the  lead  and  wing  aircraft.  In  all  channels  the  lead  responds  to  the 
formation  level  commands  while  the  wing  responds  to  the  separation  commands. 
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Dargan  iteratively  improved  the  controller  design  in  stages.  First,  he  choose  a  series 
of  test  cases  to  use  for  measuring  performance  between  each  design.  This  consisted  of  cases 
for  changes  in  the  x,  y,  and  z  separation  distances  and  cases  for  changes  in  the  velocity, 
heading,  and  altitude  of  the  formation.  Then,  he  analyzed  the  open  loop  response  and  used 
this  as  the  baseline  for  all  subsequent  designs.  The  open  loop  control  response  verified  that 
the  system  required  some  form  of  feedback.  The  formation  maneuver  commands  could  be 
followed  with  zero  steady-state  errors,  but  the  formation  separation  distances  could  not  be 
maintained. 

Dargan  then  implemented  state  feedback  for  the  formation  commands  of  velocity, 
heading,  and  altitude.  A  cascade  compensator  was  not  included  in  the  feedback  loop. 
The  longitudinal  channel  response  for  a  change  in  the  separation  distance  resulted  in  no 
steady-state  error.  The  response  for  a  change  in  velocity  resulted  in  a  steady-state  error  in 
the  X  separation  distance.  Therefore,  the  longitudinal  channel  required  some  compensation 
due  to  the  steady  state  error  from  a  change  in  velocity.  The  lateral  channel  response  was 
similar  to  the  response  of  the  longitudinal  channel.  A  change  in  y  separation  resulted  in 
no  error.  Heading  change  caused  a  steady-state  tracking  error  in  the  y  separation  distance. 
Again,  some  form  of  compensation  was  needed  for  the  lateral  channel  as  well.  The  vertical 
channel  had  a  steady-state  error  for  both  a  change  in  vertical  separation  distance  as  well 
as  a  change  in  altitude.  Compensation  was  also  required  for  this  channel.  Overall,  state 
feedback  was  not  enough  and  some  form  of  compensation  was  required  for  each  channel. 

Dargan  next  added  a  Proportional  plus  Integral  (PI)  compensator  for  control  of 
separation  errors.  This  was  composed  of  three  independent  PI  controllers;  one  for  each 
channel.  This  controller  resulted  in  a  zero  steady-state  separation  error  for  a  change  in 
the  separation  distance  and  a  change  in  formation  maneuver  for  each  channel.  However, 
the  responses  for  the  longitudinal  and  lateral  channels  had  undesirable  overshoots.  Dargan 
labeled  this  controller  as  the  unmixed  controller.  The  unmixed  controller  equations  are 

=  {Kxj,)^xe  +  (Kxi)  f  AxEdt  (2.25) 

=  {KYj,)^yE  +  (Kyi)  f  ^ysdt  (2.26) 

Jo 
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where. 


Figure  2.10  Dargan  FFCS  Block  Diagram 


hwc{i)  =  {Kzj,)Aze  +  (Kzj) 


AzEdt 


(2.27) 


•  Vwc  is  the  wing  aircraft  velocity  command  from  the  PI  controller 

•  'ipwc  is  the  wing  aircraft  heading  command  from  the  PI  controller 

•  hwc  is  the  wing  aircraft  altitude  command  from  the  PI  controller 

•  Axe^  AyE  and  Aze  are  the  respective  x,  y,  and  z  separation  errors. 

•  Kxp  and  Kxi  are  the  x  channel  proportional  and  integral  gains 

•  Kyj,  and  Ky^  are  the  y  channel  proportional  and  integral  gains 

•  Kzp  and  Kz^  are  the  wing  aircraft  z  channel  proportional  and  integral  gains. 

Lastly,  Dargan  developed  a  controller  that  used  a  mixture  of  separation  errors  and 
formation  errors.  This  was  a  mixed  controller  as  opposed  to  the  previous  unmixed  con¬ 
troller.  The  mixed  controller  included  a  linear  mixer  cascaded  with  the  PI  controller  as 
seen  in  Figure  2.10.  This  controller  worked  the  same  way  as  any  other  PI  controller.  It  at¬ 
tempts  to  null  out  an  error  between  the  commanded  input  and  the  system  response.  The 
unmixed  controller  attempted  to  null  out  the  error  between  the  commanded  separation 
distances  and  the  actual  separation  distance.  The  linear  mixer  combines  the  separation 
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distance  errors  and  the  maneuver  errors  for  each  channel  and  feeds  this  mixed  error  signal 
to  the  PI  controller  to  be  nulled  out.  The  separation  errors  for  the  longitudinal  and  lateral 
channels  are  replaced  with  the  mixed  error  signals  below 


ex  = 

(2.28) 

€y  =  k^^pE  +  kyAyE 

(2.29) 

cz  =  kzAzE 

(2.30) 

where: 

•  ex  is  the  longitudinal  channel  input  to  the  PI  controller 

•  kv  is  the  gain  on  the  velocity  error  signal,  Ve 

•  is  the  gain  on  the  x  separation  error  signal,  Axe 

•  ey  is  the  lateral  channel  input  to  the  PI  controller 

•  k^  is  the  gain  on  the  heading  error  signal,  'tpE 

•  ky  is  the  gain  on  the  y  separation  error  signal.  Ays 

•  ez  is  the  vertical  channel  input  to  the  PI  controller 

•  k^  is  the  gain  on  the  z  separation  error  signal,  Aze 
The  mixed  control  laws  are 


Vwc{i)  =  Kx^ex  +  Kx^  f  ^xdt 

Jo 

(2.31) 

i’wcii)  =  Ky^ey  +  Ky^  /  Cydt 

T  Jo 

(2.32) 

hwci'^)  —  Kzj^ez  +  Kzi  f  ezdt 

Jo 

(2.33) 

The  linear  mixer  and  PI  controller  are  shown  in  Figure  2.11.  The  x-  and  y-channels 
are  identical  except  the  Ve  term  is  replaced  with  the  term  and  Aye  replaces  Axg. 
The  only  difference  between  the  z-channel  and  the  x  and  y-channels  is  the  z-channel  does 
not  have  an  altitude  error  term  since  the  z  separation  and  the  altitude  difference  are 
synonymous. 
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UNEAR  MIXER  PI  CONTROLLER 


Figure  2.11  Dargan  FFCS  Lateral  channel 

Through  iteration,  Dargan  determined  the  gains  to  achieve  the  best  performance  [3], 
The  longitudinal  channel  still  had  an  overshoot /undershoot,  but  the  magnitude  of  the  over¬ 
shoot/undershoot  was  reduced  and  the  response  was  faster.  The  lateral  channel  resulted 
in  similar  improvements.  The  magnitude  of  the  undershoot /overshoot  was  reduced  and  a 
faster  time  response  was  achieved.  The  mixed  controller  achieved  the  desired  response  and 
concludes  the  controller  designs  developed  by  Dargan^. 

2, 6  Summary 

In  this  chapter,  the  composite  F-16  class  flight  control  system  for  the  Lead  and  the 
Wing  aircraft  was  developed.  The  large  formation  flight  controller  design  process  was 
discussed.  Finally,  the  mixed  formation  flight  controller  to  be  applied  to  this  research  was 
developed. 


^Please  note  Dargan’s  kx  and  ky  gains  are  negative  for  a  negative  feedback  system  model. 
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IIL  Upwash  and  Sidewash  Effects  on  Formation 


This  chapter  outlines  the  method  used  to  determine  the  upwash  and  sidewash  created 
by  the  Lead’s  vortex  on  the  Wing  aircraft.  The  upwash  and  sidewash  are  then  used  to 
determine  the  new  close  formation  stability  derivatives.  The  new  close  formation  stability 
derivatives  are  then  applied  to  the  Wing  FCS. 

3.1  Upwash  and  Sidewash  Derivations 

An  aircraft  moving  through  the  air  creates  vortices  behind  the  wing.  These  vortices 
can  be  seen  in  Figure  1.1  and  are  repeated  here  in  Figure  3.1.  These  vortices  exert  aerody¬ 
namic  forces  on  the  Wing  aircraft  in  the  formation.  The  formation  geometry  is  determined 
by  the  Lead’s  x,  y,  and  z  position  relative  to  the  Wing  aircraft,  in  a  rotating  frame  of 
reference  attached  to  the  Wing  as  shown  in  Figure  1.2. 


b 


Figure  3.1  Horse  Shoe  Vortex 
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Figure  3.2  Field  Strength  at  Point  P  due  to  Filament  A 

3. LI  BioUSavart.  This  section  outlines  the  method  used  to  calculate  the  upwash 
and  si  dewash  created  by  the  Lead  aircraft  on  the  Wing  aircraft.  The  derivation  of  the  effect 
of  the  Lead  aircraft’s  wing  vortex  on  the  wingman’s  wing  is  determined  by  using  an  analogy 
with  electric  field  strength  produced  by  an  electric  current  in  electromagnetics.  A  vortex 
filament  from  Fig.  3.1  is  redrawn  in  Fig.  3.2.  The  angles  /3i  and  P2  determined  by  the 
vortex  filament  segment  [A’  B’]  and  the  point  P  are  shown  in  Fig.  3.2.  The  fluid  dynamics 
analogue  of  the  Biot-Savart  law  from  electromagnetics  states  that  the  induced  velocity 
W  from  vortex  filament  A  is  given  by 

W  =  “^^“(cos/?i  -  cos/?2)  (3.1) 

47r  7*  ^ 

where  the  vortex  strength  per  unit  length  is  F,  r<;  is  the  distance  from  the  vortex  filament 
to  the  point  P,  and  the  velocity  of  the  air  wash  caused  by  the  vortex  filament  at  point  P 
and  out  of  the  page  is  W.  The  unit  vector  $  is  orthogonal  to  the  radius  vector  rc.  When 
point  A’  is  at  -oo,  the  angle  /?i  =  0,  and  when  point  B’  is  at  +oo,  the  angle  /Sj  =  tt. 
Asssuming  that  the  longitudinal  separation  in  the  formation  is  large,  i.e.,  greater  than  two 
wingspans,  is  tantamount  to  saying  that  point  B’  is  effectively  at  +oo.  When  the  actual 
distance  between  the  tail  of  the  Leader  and  the  nose  of  the  Wing  is  greater  than  2b,  where 
b  is  the  wingspan,  this  effectively  places  point  B’  at  +oo.  Then  the  Lead  aircraft’s  spanwise 
contribution  to  the  flow  field  near  the  Wing  aircraft  can  be  neglected.  Similarly,  the  Wing 
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Figure  3.3  View  from  above  the  two  aircraft 

aircraft’s  influence  on  the  flow  field  near  the  Lead  aircraft’s  wing  can  be  neglected.  These 
are  the  main  assumptions  used  in  this  research  for  the  location  of  the  Wing  aircraft’s  wing 
with  respect  to  the  vortices  caused  by  the  Lead  aircraft’s  wing.  The  vortex  goes  back  a 
great  distance  before  it  is  diffused,  validating  the  assumption  that  point  A’  is  essentially 
at  —00.  This  reduces  Eq.  3.1  to  the  simplified  form  of  the  Biot-Savart  law, 


W  = 


dr 

27rr<. 


(3.2) 


The  wash  vector,  at  the  point  P  for  the  geometry  of  Fig.  3.2,  is  out  of  the  page.  This 
would  be  the  special  case  of  no  sidewash  and  all  upwash  at  point  P  on  the  page. 

Fig.  3.3  is  a  view  from  above  the  two  aircraft  which  contains  both  the  Lead  aircraft 
replaced  by  its  horse  shoe  vortex  approximation,  with  the  sides  of  the  horse  shoe  repre¬ 
sented  by  the  A  and  B  filaments.  The  Wing  aircraft  is  represented  by  its  elliptical  wing 
approximation.  The  reduced  wingspan  of  the  Wing  aircraft  is  represented  by  V ,  This  is 
the  corrected  value  for  an  elliptical  lift  distribution  on  a  wing  and  is 


h'  =  ^-h  (3.3) 

Fig.  3.4  shows  the  view  from  behind  the  two  aircraft.  The  radius  vectors,  f,,,  from  each 
vortex  filament  to  the  wing,  disregarding  the  x  component  due  to  the  infinite  length  of  the 
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Figure  3.4  Behind  view  showing  field 


filaments,  are  for  filaments  A  and  B 


b' 

ta  =  {y-  -^-y)y -Viz-  z)z 

(3.4) 

Pb  =  (y  +  ^  -  2/)y + (^  -  z)z 

(3.5) 

where,  y  and  z  are  the  unit  directional  vectors.  Based  on  Figures  3.4  and  3.5  and  with  a 
bit  of  geometrical  insight,  it  can  be  shown  that  the  total  induced  upwash  Wuw  at  station 
(y,0)  from  both  filaments  emanating  from  the  Lead’s  wing  is 


W, 


uw 


Fi 

27r 


[(y  -  T  -  y) 


[(y  +  T  -  y)] 


\/{y  -  j  -  vY  +  y{y  +  j-yY  +  z^\_ 


<y<h' 


(3.6) 

where  x  is  the  separation  distance  in  the  x  direction,  y  is  the  separation  distance  in  the  y 
direction,  and  z  is  the  separation  distance  in  the  z  direction.  The  total  induced  sidewash 
Vsw  at  station  (0,z)  from  both  filaments  is 


V  - 


[( J  -  zf] 


y-,Q>z>-h,  (3.7) 
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3J,2  Average  Upwash  and  Sidewash.  The  average  induced  upwash  on  the  Wing 
aircraft’s  wing  is  thus  calculated  by  integrating  Eq.  (3.6): 


w  Ti 


2 

^(y  +  ^  -  y) 

V(j/  -  T  -y)"  + 

■\/(y  +  T-y)"  +  2:^J_ 

dy(-z)  (3.8) 


After  changing  the  variables  to  u  =  {{y  —  j  +  2^)  for  the  first  term  and  similarly 
w'  =  ((y  +  Y  “  vY  +  second  term,  and  changing  the  limits  to  match  with  the 

new  variables  u  and  u’,  the  integral  in  Eq.  (3.8)  is  evaluated  as 


W, 


uw. 


4x6' 


In 


{y  -  h'Y  + 


-In 


{y^b'f  +  z'^ 


+  z^ 


(3.9) 


Similarly,  averaging  by  integrating  from  the  bottom  of  the  tail  at  0  to  the  top  of  the  tail 
at  -hji  yields  the  average  sidewash,  VsWavgi  the  vertical  tail 


In 


(y  -  'if  + 

{y  -  +  (z  +  h,y 


-In 


iy+'iY  +  ^^ 

(y^  +  ^i)  +  (z  +  h,y 


(3.10) 


3.1.3  Corrected  Average  Upwash  and  Sidewash.  A  correction  term,  needs  to 
be  included  in  both  the  numerator  and  denominator  of  each  term  contained  in  the  natural 
logarithms  in  Eqs.  (3.9)  and  (3.10)  in  order  to  take  into  account  physical  viscosity  effects 
and  to  make  the  mathematical  derivation  more  accurately  follow  the  experimental  data 
[13], [14].  This  is  a  dimensionless  number  and  requires  both  the  numerator  and  denominator 
terms  inside  the  logarithm  to  be  non-dimensionalized.  To  make  the  terms  non-dimensional, 
both  the  numerators  and  denominators  are  divided  by  6^.  The  correction  term,  is  then 
added  to  both  numerators  and  denominators.  The  resulting  equations  are 


Wr 


UWa.g  - 


4x6 


In 


yl2  _|_  ^12  _j_  ^2 

{y'  -  +  iJ’^ 


-In 


(y'  +  ^r  +  z'^  +  P^ 

yl2  ^  gl2  ^  pfi 


(3.11) 
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AD 


and 

^{y'-  1)2  +  {z'  +  ^)2  +  ^2  (y'  +  1)2  +  {z'  +  '^)2  +  //2  J  ^ 

where  now  the  nondimensional  formation  geometry  parameters  are  2/'  =  f  and  z’  = 


V, 


SWava  “ 


(3.12) 


Z 


h' 


3,2  Calculation  of  Change  in  Lift  and  Drag 

The  upwash  on  the  wing  causes  a  change  in  the  angle  of  attack  of  the  wing.  This 
causes  a  rotation  in  the  lift  and  drag  vectors. 


3.2,1  Change  in  Drag.  In  Fig.  3.5,  V  is  the  velocity  of  the  aircraft,  W  is 
the  upwash,  and  V’  is  the  composite  velocity  of  the  air  at  the  surface  of  the  wing.  The 
original  lift  and  drag  vectors  are  represented  by  L  and  D  and  the  rotated  lift  and  drag  are 
represented  by  L’  and  D’,  respectively.  It  can  be  seen  from  Fig.  3.5  that  the  change  in 
angle  of  attack  is 


Aa  =  arctan 


\Wuw\ 

V 


(3.13) 


where  the  approximation  for  small  angles  has  been  applied.  It  can  also  be  seen  from  Fig. 
3.5  that  consequently  the  change  in  drag,  AD^vjdue  to  the  rotation  of  the  lift  vector  is 


ADw  =  —Lw  tan  Aa  «  —Lw 


|Waw|\ 
^  ) 


X 


(3.14) 
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where  again,  the  small  angle  approximation  is  used.  Dividing  both  sides  by  the  dynamic 
pressure  q  and  the  surface  area  S  of  the  wing,  the  non-dimensional  coefficient  of  drag 
increment  ACjy^  is 


_  ADw  _  Lw  \Ww\  _  ^  Ww 

~  qS  ~  qS  V  ~  V 


(3.15) 


Also,  the  vortex  strength  per  unit  length  can  be  expressed  as 


Ll  _  Xi 

pviy  pv\b  pv\h 


'kA, 


-Ci.Vb 


(3.16) 


where  S  =  area  of  wing,  b  =  wing  span,  V  =  airspeed  of  formation,  =  aspect  ratio  of 
wing,  and  Cl^  =  lift  coefficient  of  the  Lead  aircraft.  Substituting  the  above  equation  for 
vortex  strength  and  also  substituting  the  previously  derived  upwash  expression  into  the 
change  in  drag  equation,  yields 


ACn.„  = 


r  A 


+  +  p'' 


In 


(2/'  -  |)^  +  z'"^  A  pi^ 


-In 


[y' + 

yl2  ^  ^  y2 


(3.17) 


3.2.2  Change  in  Lift.  The  change  in  lift,  is  given  by 


ACij,,  —  A.aaw  = 


\Wuw\ 

V 


o-w 


(3.18) 


where  aw  is  the  lift  curve  slope  of  the  wing.  Substitution  of  the  upwash  and  vortex  strength 
values  as  before,  results  in 


ACi,,= 


^  2 
jAr  TT^ 


In 


yl2  _|_  ^/2  _|.  y2 

{y'  -  f)^  +  z'^  +  p^ 


-In 


{y' A  \f  +  Z'^  +  p-^ 


yl2  ^  ^12  ^  ^2 


(3.19) 


This  is  the  change  in  lift  coefficient  of  the  Wing  from  the  upwash  created  by  the  Lead 
aircraft. 
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3,2,3  Change  in  Side  Force,  The  sidewash  created  by  the  Lead  also  causes  a 
change  in  the  force  on  the  vertical  tail.  This  change  in  side  force,  Ai^y,  is 


AFy  = 


(3.20) 


where  ri  is  the  aerodynamic  efficiency  factor  at  the  tail,  S^t  is  the  vertical  tail  area,  and 
a^t  is  the  lift  curve  slope  of  the  vertical  tail.  The  non-dimensional  change  in  side  force 
coefficient,  is 

(3.21) 


Ari  r5w| 

ACy  =  V—ayi  ^ 


Inserting  the  sidewash  expression  of  Eq.  (3.12)  back  into  the  non-dimensional  side  force 
Eq.  (3.21)  gives  the  side  force  coefficient 


ACy  =  7] 


Syt  O-vt  Tx 

S  V  4^ 


(3.22) 


In 


(y'  -  1)^  + 


-In 


/+|)2  +  ^/2+^2 


{y'  -  1)2  +  (^/  -P  1^)2  +  fjfl  (y>  +  1)2  +  (^/  +  1^)2  +  ^2  J 


Finally,  substituting  back  in  for  F  by  using  Eq.  (3.16),  the  change  in  side  force  coefficient 
becomes 


ACy  = 


1  r]SytaytCL^b2 


ttAi 


2SK 


(3.23) 


In 


{y'  - 1)'  +  2'"  + 


In 


{y'  + 1)'  +  +  m" 


{y'  -  1)2  +  {z'  -f  ^)2  +  ;,2  +  1)2  +  +  i^)2  +  ^2  J 
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3.3  Modified  Wing  Aircraft  Control  System 

It  is  envisaged  that  each  aircraft  is  equipped  with  a  flight  control  system  that  includes 
three  standard  autopilots:  Heading  hold,  Mach  hold,  and  Altitude  hold  autopilots. The 
aircraft  autopilots  that  makeup  the  FCS  are  developed  in  Chapter  II.  The  formation  flight 
control  autopilot  resides  on  the  Wing  aircraft.  It  is  an  outer-loop  controller  which  receives 
measurements  of  the  Lead  aircraft’s  position  relative  to  the  Wing  aircraft  and  it  drives 
the  reference  signals  of  the  Wing’s  three  axes  Mach  hold.  Heading  hold,  and  Altitude  hold 
autopilots. 

The  formation  flight  control  autopilot  resides  on  the  Wing  aircraft.  It  is  an  outer-loop 
controller  which  receives  measurements  of  the  Lead  aircraft’s  position  relative  to  the  Wing 
aircraft  and  it  drives  the  reference  signals  of  the  Wing’s  three  axes,  Mach  hold.  Heading 
hold,  and  Altitude  hold  autopilots. 

For  close  formation  flight,  the  Wing  FCS  needs  to  be  modified  to  account  for  the 
additional  aerodynamic  interactions  created  by  the  upwash  and  sidewash  from  the  Lead 
aircraft.  The  upwash  causes  an  aerodynamic  force  in  the  x  and  z  directions  in  the  form  of 
a  change  in  drag  and  change  in  lift  as  derived  above.  Thus,  the  Wing  aircraft  needs  to  be 
retrimmed  in  pitch.  The  sidewash  induces  a  force  in  the  y  direction  caused  by  a  change 
in  lift  on  the  vertical  tail,  which  requires  the  Wing  aircraft’s  lateral  directional  control 
channel  to  be  retrimmed.  This  change  in  forces  need  to  be  trimmed  out  by  the  Wing’s 
FCS.  Above  and  beyond  this  retrimming  action,  it  is  important  to  include  the  dynamical 
change  in  forces  caused  by  a  perturbation  Ax,  Ay,  and  Az  in  the  Lead’s  position  in  the 
formation  relative  to  the  Wing. 

The  upwash  and  sidewash  also  cause  changes  in  moments  applied  to  the  wing  aircraft. 
However,  moment  effects  are  not  explored  in  this  paper  because  of  the  rather  rudimentary 
modeling  employed,  where  the  Wing  aircraft  is  considered  a  point  mass. 

Thus,  new  stability  derivatives  for  change  in  Wing  drag  resulting  from  a  change  in 
X,  y,  and  z  positions  in  the  formation  need  to  be  determined  for  the  Mach  hold  channel. 
Similarly,  the  new  stability  derivatives  for  the  change  in  lift  due  to  a  change  in  the  x,  y, 
and  z  positions  need  to  be  determined  for  the  Altitude  hold  channel.  Finally,  new  stability 
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derivatives  for  the  change  in  side  force  due  to  a  change  in  x,  y,  and  z  positions  in  the 
formation  need  to  be  determined  for  use  in  the  Heading  hold  channel.  The  resulting  Wing 
FCSis 


/yy 


hyy 


It.  It. 

- TV  H - TVc  + 

Tv  Tv 

^  +  ACoyy^y  +  ACd^,z] 

1  /  1  , 

- Ww  d - i^Wc  + 

T~%fj  T'tj} 

+  ^Cy^^y  +  ACy^^z] 


m 


[ACl^^x  +  AC^^^y  +  ACl^^z 


where:  V\y^  is  the  reference  signal  to  the  Wing  Mach-hold  autopilot,  and  similarly,  hw^  is 
the  reference  signal  to  the  Altitude-hold  autopilot  and  is  the  reference  signal  to  the 
Heading- hold  autopilot.  Vw  is  the  Wing’s  velocity,  'i/^w  is  the  Wing’s  heading,  hw  is  the 
Wing’s  altitude,  and  and  z  are  the  perturbations  in  the  Lead’s  position  in  the  formation 
relative  to  the  Wing  from  the  nominal  location  (x^y^z).  In  the  V\y  differential  equation 
the  new  stability  derivatives  for  the  longitudinal  and  the  vertical  perturbations,  ACjd^^, 
and  ACd^^^  and  ACl^^^  and  ACl^^  are  multiplied  by  Multiplying 

them  by  qS  converts  them  back  to  a  force,  and  then  dividing  by  the  mass  converts  them 
to  an  acceleration,  as  required.  The  new  stability  derivatives  for  the  heading  differential 
equation  are  multiplied  by  The  additional  division  by  the  velocity  is  required  because 
the  heading  rate  is  an  angular  velocity.  Also,  the  centrifugal  acceleration  is  A  =  u;  X  V 
where  V  is  the  aircraft’s  velocity  and  u;  is  the  angular  velocity.  Since  all  of  these  vectors 
are  orthogonal  to  each  other,  u  =  AjV .  Hence,  the  new  stability  derivatives  in  the  heading 
equation,  when  multiplied  by  ^  are  converted  to  an  angular  velocity,  as  required  in  the 
heading  equation. 


3,3,1  Calculation  of  the  New  Formation  Stability  Derivatives,  The  change  in 
drag,  lift,  and  side  force  on  the  Wing  aircraft  in  a  close  formation  have  been  previously 
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calculated.  To  determine  the  change  in  these  forces  due  to  a  change  in  the  x,  y,  and  z 
Lead’s  relative  position  in  the  formation,  a  linearization  is  performed  about  the  nominal 
Lead  position  (which  is  measured  with  respect  to  the  Wing’s  position)  in  the  optimal  close 
formation:  The  latter  is  y  =  and  z  =  0.  This  requires  derivatives  of  the  change  in  drag, 
lift,  and  side  force  to  be  evaluated  at  these  values  for  y  and  z. 

First,  a  dimensionless  expression  is  introduced  for  the  change  in  induced  drag,  (Juwiy'i  ^')'- 


2  (y/  +  1)2  +  ^/2  +  ^2 

’T'  r  (y'  -  V'  +  y'^  + 


(3.24) 


A  similar  non-dimensional  function  is  defined  for  the  sidewash  component, 


(Tswiy',  z')  = 


In. 


{y' -  ly  +  z'^  + 


—  In  • 


{y' +  ly  +  z'^  + 


(y'  - 1)^  +  +  ^y  +  (y'  + 1)^  +  (^'  +  ^y  + 

Based  on  these  definitions,  the  change  in  lift,  drag,  and  side  force  are  expressed  as 


(3.25) 


1 


ClL^Lw'^Uwiy'  1  z') 
O’CLw^Uwiy'iZ') 

CiL^swiy'  ,z') 


T^Ar 
t^Ar  2Shg 


(3.26) 

(3.27) 

(3.28) 


The  only  parts  of  the  above  equations  that  vary  with  a  change  in  the  x,  y,  and  z  position 
are  the  (Juwiv' iz')  function  for  the  change  in  drag  and  lift,  and  the  <Tsw(y',z')  function 
for  the  change  in  side  force.  The  derivatives  of  cTuw{y'yz')  and  (Tsw{y',z')  are  analytically 
obtained  as  follows:  The  partial  derivatives  for  (Tuw(,y'yZ')  are  calculated  as 


d(7\ 


UW  I 


dx'  ' 
duuw  I 


da. 


UW  I 


dz' 


o 

II 

o 

Jl 

N 

II 

(3.29) 

(3/8)7r 

=  [(£)2  +  ^2][(^)2  +  ^2] 

(3.30) 

li 

\\ 

o 

II 

O 

(3.31) 
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and  the  partial  derivatives  for  (Tsw{y',z')  are  calculated  as 


da, 


5VV  I 


dx' 


ly' 


=  f,^'=o  =  0 


daswf  _ 

dy' 


l(W 


i(fF+»‘’][(i)’+(^)»+/.’i 

3/^\2 


2\  b  ) 


da 


siv 


dz' 


|y'=f,^'=0  - 


[(|)2+^2  +  (^)2]  f(4)2  +  ^2  +  (^)2] 


(3.32) 

(3.33) 

(3.34) 


Inserting  the  above  derivatives  into  the  change  in  drag,  lift,  and  side  force  coefficients 
expressions  yields 


AC. 


Dwx 


Lwx  ^^YvVa 

1 


-c 


1 


AO 


Ywy  “ 


1 


TrAf 


4Sb 


-  ACDw,  =  ACLyy,^  =  0 

(3/8)7r 

(3.35) 

(3/8)lr 

(3.36) 

(3.37) 

l(f )’  +  «“]  [(«’  +  (t^)’  +  «’l  [(f)"  +  M"1  [(¥)"  +  (f )"  +  ^"j . 


\  1  V^vt^vt  ^ 

^^YWz  = - ;;; - rr; — 

ttAr  25 


(fKi) 


(3.38) 


3.3.2  Modified  Wing  Aircraft  Control  System.  Based  on  the  above  derived  sta¬ 
bility  derivatives  associated  with  the  forces  created  by  the  upwash  and  sidewash,  the  new 
Wing  FCS  is 


% 

hw 

i’w 


~Vw  +  f  W,  +  §ACd„,!, 

Ty  Ty 


m 

1  1  . , 

H - I  f^w  — 

/  ^ha^hij 

1  ,  .  1 


=  -(—  +  —) 
\rh,  ThJ 


1  L  1  ..  qS 

hw  +  ——hwo  +  —ACL^y 
naT^h,,  m 


'tl^w  +  - — V’We  +  ^  [ACy^^^y  +  ACy^,z] 


T^/iv 


mV 


(3.39) 

(3.40) 

(3.41) 
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where  ACd^^,  ACyw^,  and  ACy^,  are  the  new  close  formation  stability  deriva¬ 

tives. 

3.4  Summary 

The  average  upwash  and  si  dewash  effects  from  the  Lead  vortex  on  the  Wing  aircraft’s 
wing  were  developed.  These  effects  were  then  used  to  derive  the  change  in  drag,  change  in 
lift,  and  change  in  sideforce,  AD,  AX,  and  AT,  respectively.  New  close  formation  stability 
derivatives  were  developed  and  the  Wing  FCS  was  modified  to  account  for  the  additional 
aerodynamic  forces  exerted  by  the  Lead’s  vortex  on  the  Wing. 
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IV.  Control  Design  and  Evaluation 


The  complete  linear  and  nonlinear  system  models  are  developed  in  this  chapter.  The 
method  for  determining  the  best  gains  for  the  close  FFC  is  described.  Finally,  a  stability 
and  robustness  analysis  is  performed  for  the  complete  system  using  these  gains. 


4.1  Complete  System  Model 

Equations  (3.39)  -  (3.41)  are  adjoined  to  the  nonlinear  kinematics  Equations  (2.22)  - 
(2.24)  where  Equation  (3.41)  for  ^  is  inserted  into  the  x  and  y  separation  differential  equa¬ 
tions,  yielding  the  hybrid  nonlinear  six-dimensional  close  formation  flight  control  system: 


X  = 


y  = 


Vw 

i’w 


y  y 

- -0^  --  _l_  Vj^cos('ipL  —  '^w)  H - + 

'^tpw  '^'^w 

X  X 

- i>w  +  VLsin{ipL  -i^w) - V’w.  - 

-—Vw  +  —Vw^  +  -  y) 

Ty  Ty  TTl 

-- — +  - — tAwc  +  [ACrwyiy  -y)  +  -  ^)] 

TnVw 


z  =  C 


;  1  1.  1, 

C  —  —  ( - "h  —  I  C  — - ^  T - l^w  — - l^L  "h 

VT-ft.  ThJ  Th,Th^  Th^Thy  °  Th^Thy  ° 


^AC,^^{y-y) 


(4.1) 


(4.2) 


(4.3) 

(4.4) 

(4.5) 

(4.6) 


where  the  vertical  separation  2:  =  hw  —  and  the  Wing  and  Lead  aircraft  are  assumed 
to  have  the  same  vertical  dynamics. 

The  six  states  are  x^y^'ipw^Vw^z^  and  The  three  control  inputs  to  the  respective 
Wing’s  Heading  hold,  Mach  hold,  and  Altitude  hold  autopilot  channels  are  'il)wc^  Vvci 
and  hwc-  The  Lead’s  control  inputs  are  viewed  as  a  disturbance;  thus,  the  disturbance 
signals  are  and  hi^.  This  is  the  full  nonlinear  model  used  in  the  simulation.  For 

the  purpose  of  controller  design,  linearization  of  the  above  hybrid  system  yields  the  linear 
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perturbation  equations 


-Vw  +  Vl-^ — —i^Wc  +  y  +  ACy^  z  (4.7) 

mV  I  ^  *  J 

-v\ij)w  ^V'lj^L-  —i^Wc  “  WcY^  y  +  ACy  ^:]  (4.8) 

+  +  (4.9) 


7^^*^  +  '*’  mF  (4-10) 

(4.11) 

^  ~  -  :r4“^^e  (4-i2) 


'^ha  '^hh 


'^htt'^hb  '^ha'^hb 


The  resulting  state  space  representation  based  on  the  above  defined  states,  controls, 
and  disturbances  is 


d  y 

dt  j. 


+  B  '4’wc  +  r  V’L 


(4.13) 


where  the  dynamics  matrix 
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^ACy^^y 

^ACdw^ 


■t^g 


0  0 
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0  0 
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■^^^YwyX  (7^“^i)G  --^ACy^^X 

^ACy^^llG  :^ACy^^llG 


(^  +  i)  J 

(4.14) 
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the  input  matrix 


0 

-^G 

0 

1 

TVyy 

0 

0 

0 

X  Q 

'^i>w 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

and  the  disturbance  matrix  is 


10  0 


r  = 


0  0  0 

0  VlG  0 

0  0  0 

0  0  0 


(4.15) 


(4.16) 


and  where  the  number  G  =  tt/ISO  and  is  used  to  keep  consistency  among  units.  The  block 
diagram  for  the  above  linearized  system,  Equations  (4.14)  -  (4.16)  is  shown  in  Figure  4.1. 


The  dynamics  A  matrix  shows  that  there  is  additional  coupling  introduced  to  the 
system  above  and  beyond  the  kinematically  induced  coupling.  This  is  especially  evident 
in  the  term  which  couples  the  y  separation  into  the  altitude  hold  autopilot  and  the 
A3_5  term  which  couples  the  z  separation  into  the  heading  hold  autopilot.  The  heading 
is  kinematically  coupled  into  the  y  separation  equation  and  the  y  separation  is  aerody- 
namically  coupled  into  the  Mach-hold  autopilot,  terms  A3_5  and  Ai_3  respectively.  This 
leads  to  an  indirect  coupling  between  the  z  separation  and  the  Mach-hold  autopilot.  There 
is  no  coupling  between  the  horizontal  and  vertical  channels  for  a  large  formation.  This 
interaction  between  channels  is  one  of  the  most  interesting  effects  of  close  formation  flight 
and  is  discussed  further  in  Chapter  V. 
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Figure  4.1  Linearized  System  and  Formation  Flight  Controller 

Typical  characteristics  for  an  F-16  class  aircraft  are  listed  in  Table  4.1  for  an  altitude 
of  45,000  ft,  velocity  of  0.85  Mach,  and  a  dynamic  pressure  of  155.8  Ibsfft^.  The  close 
formation  stability  derivatives  for  this  flight  condition  and  the  associated  characteristics 
are  listed  in  Table  4.2. 

The  dynamic  A  matrix  evaluated  for  the  typical  characteristics  listed  in  Table  4.1 
and  the  nominal  x,  y,  and  z  separations  listed  above,  with  no  close  formation  induced 
aerodynamic  coupling  effects,  is 


0  -1.0000  0  -0.3084  0  0 

0  -6.0000  0  0  0  0 

0  0  0  -13.6136  0  0 

A  =  (4.17) 

0  0  0  -0.7500  0  0 


0  0  0  0  0  1.0000 

0  0  0  0  -0.8447  -3.5118 
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The  dynamic  A  matrix  evaluated  for  the  data  listed  in  Table  4.1  and  the  nominal  x,  y, 
and  z  separations  with  the  close  formation  coupling  eflfects  included  is 


0 

-1.0000 

0.0057 

-0.3084 

-0.0020 

0 

0 

-6.0000 

-0.0471 

0 

0 

0 

0 

0 

-0.0144 

-13.6136 

0.0050 

0 

0 

0 

0.0138 

-0.7500 

0.0138 

0 

0 

0 

0 

0 

0 

1.0000 

0 

0 

0.4663 

0 

-0.8447 

-3.5118 

It  can  be  seen  from  a  comparison  of  Eqs.  (4.17)  and  (4.18)  that  the  coupling  terms 
introduced  by  the  close  formation  are  small.  One  notable  exception  is  Ae.a  =  0.4663.  This 
justifies  the  use  of  the  large  formation  controller  since  this  is  essentially  a  three  channel 
system,  although  special  care  is  required  in  the  z-channel  controller  design. 
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Table  4.1  F-16  Class  Aircraft  Characteristic  Values  This  data  corresponds  to  the  fol¬ 

lowing  flight  condition:  Altitude  of  45,000  ft,  dynamic  pressure  155.8  Iblft^ 


Wing  Area 

S 

300  ff 

Wing  Span 

b 

30  ft 

Aspect  Ratio 

Ar 

3 

Lift  Curve  Slope 

a 

5.3  per  rad 

Tail  Area 

s.t 

54.75 

Tail  Height 

K 

120  in 

Tail  Lift  Curve  Slope 

dvt 

5.3  per  rad 

Aerodynamic  Efficiency  Factor 

n 

0.95 

Velocity  Time  Constant 

TVw 

5  sec 

Heading  Time  Constant 

0.75  sec 

Altitude  Time  Constant 

h/Q 

0.3075  sec 

Altitude  Time  Constant 

hi 

3.85  sec 

Gross  Mass 

m 

776.4  slugs 

Gross  Weight 

w 

25000  lb 

Velocity 

V 

Mach  =  0.85  or  825  ftjsec 

Table  4.2  Close  Formation  Stability  Derivatives  for  y  =  and  z  =  0 


-0.000782 

-0.0077 

AC'y,,„ 

0.0033 

-0.0011 
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4-2  Control  Design 

Similar  to  the  large  FFC  system,  note  that  the  y  and  z  channel  states  y,  ipw,  z,  and  C, 
are  decoupled  from  the  x  channel  states,  x,  and  Vy/\  which  is  easily  seen  in  Equation  4.14. 
Hence,  one  first  designs  a  controller  for  the  y  and  z  channels,  where  the  control  signals 
are  V’Wc  following  which  one  turns  to  the  design  of  the  x-channel  controller,  viz., 

one  synthesizes  the  control  law  for  the  control  signal.  Furthermore,  even  though  the 
(y,z)  model  does  not  further  decompose  into  individual  y  and  z  channels  as  in  [8]  where 
the  coupling  is  exclusively  induced  by  the  kinematics  of  the  large  FFC,  the  same  controller 
will  be  applied  because  the  close  formation  induced  coupling  is  weak. 

The  controller  contains  a  linear  mixer  on  the  x  and  y  channels  and  a  standard  PI 
compensator  on  the  x-,  y-  and  z-channels.  The  FFC  Equations  2.31  -  2.33  developed  in 
Chapter  II  are  repeated  below. 


Vwc(t)  = 

tpwcit)  = 

hwcit)  = 


Kx^lkvVs  +  k^AxE]  +  Kxj  J  [kyVE  +  k^AxEldt 

Ky^iki^i^E  +  kyAys]  +  Ky,  [k^tpE  +  kyAyE]dt 

Kz^k^AzE  +  Kzi  f  k;,AzEdt 

Jo 


The  FFC  gains  have  been  choosen  based  on  the  three  criteria  listed  below: 


•  The  maximum  positive  and  negative  perturbations  from  nominal  for  a  single  channel 
should  be  approximately  equal  for  any  heading  maneuver  within  the  test  envelope 
defined;  e.g.,  the  y  separation  for  a  -30°  heading  change  ranges  from  a  minimum  of 
20.5  ft  to  a  maximum  of  26.3  ft  or  -3  ft  and  2.8  ft  from  nominal  respectively.  The 
deviation  to  the  right  is  approximately  equal  to  the  deviation  to  the  left. 

•  The  responses  should  be  consistent  for  either  a  positive  or  a  negative  change  in  either 
the  heading,  the  velocity,  or  the  altitude;  e.g.,  a  -30  degree  heading  change  will  cause 
a  maximum  error  of  4  ft  in  x  and  3  ft  in  the  y  and  a  -t-30  degree  heading  change  will 
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cause  a  maximum  error  of  3.3  ft  in  the  x  and  3  ft  in  the  y.  Both  maneuvers  result  in 
errors  of  approximately  the  same  magnitude. 

•  The  Wing  aircraft  should  not  cross  over  the  x  ox  y  axes.  This  criterion  has  the  dual 
effect  of  preventing  collisions  between  the  Lead  and  the  Wing;  also  in  a  formation 
containing  both  a  left  and  right  Wing  aircraft. 

Previous  AFIT  students  [6], [12]  used  SISO  techniques  to  determine  the  gains.  The 
method  consisted  of  successively  closing  loops  to  create  a  new  A  matrix.  After  each  loop 
closure,  the  new  A  plant  matrix  is  used  for  the  next  loop  closure.  Once  the  appropriate 
loops  were  closed,  the  root  locus  was  used  to  iteratively  determine  the  best  gains  to  achieve 
the  desired  performance.  Even  though  this  method  is  not  completely  rigorous  for  a  MIMO 
system,  it  does  give  a  good  indication  of  gain  magnitudes  and  system  response. 

4 -2 A  Gain  Selection  Method.  A  simple  systematic  approach  was  used  to  itera¬ 
tively  determine  the  best  gains.  First,  only  the  completely  linearized  system  as  shown  in 
Equations  4.7  -  4.12  and  in  Figure  4.1  is  considered.  With  the  cross  coupling  terms  zeroed 
out,  Dargan’s  gains  were  chosen  as  the  starting  point  [3].  The  system  is  subjected  to  a 
positive  and  negative  velocity  change  and  the  horizontal  channel  gains  are  adjusted.  The 
system  is  then  subject  to  a  positive  and  negative  heading  change  and  the  horizontal  gains 
are  adjusted.  Lastly,  the  system  is  subjected  to  a  positive  and  negative  altitude  change 
and  the  vertical  channel  gains  are  adjusted.  The  order  of  gain  adjustment  can  affect  the 
number  of  iterations  required  to  determine  the  best  gains.  The  best  order  of  gain  adjust¬ 
ment,  i.e.,  the  method  which  results  in  the  fewest  iterations,  was  to  adjust  the  mixer  gains 
first,  to  achieve  the  best  maneuver  tracking  response,  and  then  adjust  the  PI  compensator 
gains  to  achieve  the  best  formation  response.  Once  the  best  gains  are  determined  for  the 
complete  linear  system  with  no  coupling,  the  next  step  is  to  add  back  the  linear  cross 
coupling  terms.  Finally,  apply  these  gains  to  the  full  non-linear  model  shown  in  Figure 
4.17  and  adjust  the  gains  for  the  best  response. 

4^2.2  Gain  Selection.  The  Six  cases  used  for  gain  selection  are:  +/-50  ft/sec 
velocity  change,  +/-  30°  heading  change,  and  -h/-1000  ft  change  in  altitude.  Figures  4.2 
and  4.3  show  the  responses  of  the  system  using  the  large  FFC  with  Dargan’s  gains  for  a 
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50  ft/sec  change  in  velocity.  Figures  4.4  and  4.5  show  the  system  responses  for  a  “30® 
heading  change  using  Dargan’s  gains.  Figures  4.6  -  4.7  show  the  response  for  a  1000  ft 
increase  in  altitude. 

The  format  for  the  time  response  plots  and  the  Lissajous  plots  is  described  here. 
The  time  plots  contain  six  subplots.  In  the  top  three  subplots,  the  solid  line  represents 
the  Lead  aircraft  and  the  dashed  line  represents  the  Wing  aircraft.  In  some  cases  only  a 
solid  line  is  visible,  this  is  due  to  the  overlap  of  the  responses  from  the  Lead  and  the  Wing. 
In  the  bottom  three  subplots,  the  solid  line  represents  the  nominal  separation  values  and 
the  dashed  line  represents  the  actual  values.  Starting  points  in  the  Lissajous  plots  are 
represented  by  an  ‘O’  and  the  ending  points  are  represented  by  an  ‘X’.  This  format  is  used 
for  each  plot  listed  in  the  following  sections. 

As  expected,  Dargan’s  large  FFCS  gains  resulted  in  poor  performance  for  a  close 
formation.  A  velocity  decrease  of  50  ft /sec,  Figures  4.2  and  4.3  resulted  in  a  long  delay 
before  the  Wing’s  velocity  responded  to  the  Lead’s  velocity.  This  delay  causes  the  Wing 
to  overtake  the  Lead,  but  since  the  Lead’s  wing  and  the  Wing’s  wing  overlap  and  both 
aircraft  are  at  the  same  altitude,  a  collision  occurs.  This  is  easily  seen  when  the  separation 
in  the  x  direction  decreases  to  zero,  crosses  the  y  axis,  and  becomes  negative.  A  —30® 
heading  change  resulted  in  a  maximum  deviation  of  4  ft  |b)  in  the  x  direction  and 
a  maximum  deviation  of  41  ft  («  1-  b)  in  the  y  direction,  see  Figures  4.4  and  4.5.  The 
deviation  in  the  y  direction  crossed  the  x  axis.  This  would  most  likely  have  resulted  in  a 
collision  for  the  case  of  two  Wingmen.  The  responses  to  an  increase  in  altitude  of  1000  ft 
using  Dargan’s  large  FFCS  gains  is  shown  in  Figures  4.6  -  4.7.  The  altitude  also  suffered 
from  a  long  delay  before  the  Wing  responded  to  the  Lead’s  ascent.  The  separation  in  the 
z  direction  grew  to  over  400  ft  («  13b).  No  collision  would  occur,  but  the  enormous  error 
is  still  not  acceptable  for  close  formation  flight.  Therefore,  all  three  maneuver  responses 
are  unacceptable  for  close  formation  flight  and  require  the  gains  to  be  adjusted. 
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Dargan  Gains 


Gains 


With  the  system  responses  from  Dargan’s  gains  as  the  baseline,  adjustment  of  the 
X“Channel  was  chosen  first  for  the  close  FFC.  Since,  it  is  independent  of  the  other  channels 
when  only  a  velocity  change  occurs,  it  should  be  the  easiest  of  the  two  horizontal  channels 
to  determine  the  gains.  Using  the  Dargan  gain  response  as  the  baseline,  the  system  was 
subjected  to  a  50  ft /sec  change  in  velocity.  The  velocity  error  gain  (ky)  was  increased 
to  get  the  Wing  velocity  to  track  closely  with  the  Lead.  Once,  the  best  velocity  response 
was  achieved  the  x  separation  error  gain  (k^;)  was  increased  to  get  the  best  x  separation 
response.  Both  of  these  gains  reached  a  certain  limit  in  which  further  increase  resulted 
in  degraded  performance.  Once  this  point  was  achieved  the  x  channel  proportional  gain 
(Kxp)  was  adjusted  to  give  the  best  x  separation  response.  This  gain  also  reached  a  limit 
in  which  further  increase  resulted  in  degraded  performance.  The  x  separation  looked  like 
an  elongated  bell  sitting  on  the  line  marking  the  nominal  x  value.  This  shape  indicates 
that  the  gain  on  the  integrator  was  too  small  because  the  error  continued  to  build  up  for 
a  long  time  and  then  took  a  long  time  to  decay.  Therefore,  the  x  channel  integration  gain 
(Kxj)  was  increased.  The  combined  effect  of  all  gain  adjustments  resulted  in  an  initial 
ramp  like  increase  in  x  separation  error,  but  with  a  much  smaller  magnitude,  then  a  quick 
dropoff  back  to  the  nominal  x  separation  distance.  This  was  first  done  for  an  increase 
in  velocity.  The  gains  were  then  readjusted  to  get  the  best  average  response  for  both  an 
increase  and  decrease  in  velocity  to  fulfill  the  second  gain  criterion. 

These  new  gains  provided  the  starting  point  for  the  next  phase.  This  phase  consisted 
of  inputting  a  heading  change.  This  was  a  little  trickier  since  heading  changes  affect  both 
the  y-  and  x-channels  simultaneously.  However,  applying  the  same  approach  as  above 
worked  well.  First,  the  heading  error  gain  (k^)  was  adjusted  to  track  the  Lead’s  heading 
as  closely  as  possible.  The  velocity  and  x  separation  responses  degraded  slightly  from  the 
heading  gain  adjustment  and  required  a  slight  re- adjustment  of  the  velocity  and  the  x 
error  gains  (ky  and  k^o).  The  y  error  gain  (ky)  was  then  adjusted.  Then  the  proportional 
y  channel  gain  (Kyp)  and  finally  the  integral  y  channel  gain  (Ky^)  was  adjusted.  This 
method  worked  well  and  was  rather  logical  with  respect  to  the  system  responses  for  each 
step. 
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Time  (sec) 

Figure  4.6  1000  ft  Altitude  Increase  Plot  Using  Dargan  Gains 

The  adjustments  of  the  z  channel  gains  were  almost  trivial  with  respect  to  the  x- 
and  y-channels.  It  required  the  adjustment  of  three  gains,  even  though  adjustment  of  the 
linear  mixer  gain  (kjg)  is  the  same  as  adding  an  equal  gain  to  both  the  proportional  (Kzp) 
and  integral  (Kzj)  z  error  gains.  First,  the  z  error  gain  (k^)  was  increased  to  achieve  a 
good  response.  The  proportional  (Kzp)  gain  was  then  increased  to  improve  the  response. 
The  result  of  these  two  adjustments  was  the  best  response  and  did  not  require  any  further 
adjustment  of  the  integral  gain  {Kz^)-  SISO  root  locus  methods  would  also  have  worked 
well  for  determining  the  best  gains  for  this  channel. 

The  simulations  were  then  run  again  using  the  nonlinear  kinematics  Equations  (4.1) 
-  (4.7)  and  no  cross  coupling  to  determine  if  the  gains  determined  in  the  previous  run 
needed  further  adjustments.  Adjustments  were  not  needed.  Finally,  cross  coupling  terms 
were  included;  again,  gain  adjustments  were  not  necessary.  The  final  gain  values  are  listed 
in  Table  4.3. 
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4.3  Stability 


In  order  to  determine  the  stability  of  the  closed  loop  system,  the  linearized  plant  and 
controller  must  be  combined  together  in  an  augmented  closed  loop  plant.  The  eigenvalues 
of  this  new  system  should  indicate  whether  or  not  the  complete  linear  system  is  stable.  The 
actual  formation  flight  controller  is  composed  of  nonlinear  saturation  limits  associated  with 
the  FCS.  For  this  reason,  stability  cannot  be  guaranteed,  but  for  the  operating  envelope 
deflned,  it  is  a  good  indication  of  system  stability. 

Since  the  FFCS  consists  of  a  linear  mixer  feeding  into  a  three  channel  PI  controller, 
three  more  states  are  needed  for  the  augmented  matrix.  The  three  augmented  states  are 


kv  {Vl  -  Vw)  +  k^  {xg  -  x)  =  ex 

(4.19) 

k,p  (V’L  -  M  +  ky  {yg  -y)  =  eY 

(4.20) 

II 

1 

(4.21) 

where,  a,  (3,  and  7  are  not  related  to  the  angle  of  attack,  the  sideslip  angle,  or  the  flight 
path  angle.  The  new  states  are  the  outputs  of  the  respective  integrator  for  each  channel’s 
PI  controller  as  shown  in  Figure  2.11.  The  input  commands,  Equations  (2.31)  -  (2.33)  are 
now  solved  using  these  new  states.  The  resulting  equations  are 

Vwc  —  Kxjd  +  Kxpkv  (Vl  —  ^w)  +  Kp^k,.  —  x)  (4.22) 

i>Wc  =  Ky:I3  +  Krpk^  {■^l  -  i’w)  +  Kypky  {y^  -  y)  (4.23) 

hwc  =  Kz,'f  +  Kzpk^  (zc  -  z)  (4.24) 

Equations  (4.19)  -  (4.21)  are  augmented  to  Equations  (4.7)  -  (4.12)  and  the  new 
input  Equations  (4.22)  -  (4.24)  are  substituted  into  Equations  (4.7)  -  (4.12).  This  creates 
a  new  nine  state  system.  The  new  inputs  are  Xc,  yc,  and  Zg  and  the  Lead  is  again  treated 
as  a  disturbance.  The  resulting  state  space  representation  based  on  the  augmented  states. 
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controls,  and  disturbances  is 
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where  the  is  arbitrarily  broken  into  and  Aci2,  and  Bd,  and  Tci  matrices  are, 
respectively, 

Aei  =  [AciiAd^]  (4.26) 


0 

-1 

- 

^(i  +  hry^ir,^)G  1 

Tv 

-^(l  +  iTxpi^v) 

^ACdwv 

0 

0 

0 

SKYpK^  \  ^ 

0 

0 

(  ^  +  SACy,,,(l/G))  - 

KypKy) 

0 

0 

0 

0 

0 

0 

^^Clwy 

0 

“  Ajv 

0 

0 

0 

0 

-ky 

-k 

* 

0 

0 

0 

0 

(4.27) 

0 

0 

^yG 

0 

0 

0 

^Xr 

Ty 

0 

0 

-X-^ACy^, 

0 

0 

0 

■^ACYyi,,(l/G) 

0 

0 

£il 

^V'W' 

0 

Ac12  = 

0 

1 

0 

0 

0 

(4.28) 

/  1  . 

0 

0 

Kzjk^ 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-k. 

0 

0 

0 

0 

4-16 


0 

Kxpk^ 

Tv 

0 

0 

0 

0 

0 

0 


—  Kyp  ^  y 

0 

-x^^p^y 

T^W 

Ky^Jz^ 

‘^•*w 

0 

0 

0 


0 


0 

0 

0 

0 

0 

Kzp 
’■ka  ’■*6 

0 

0 

h 


Tci  = 


1 

KxpK, 

TV 

0 

0 

0 

0 

ky 

0 

0 


_yi[lEhLQ 

^  ’’V-vv 


0 


fy_^Kv^\  Q 
\  ^V'W  / 

K  Yp  k 
T^W 


0 

0 


0 

k^ 

0 


0 

0 

0 


0 

0 


’■k.  ’■ki, 
0 

0 

0 


where  the  number  G  =  tt/ISO  and  is  used  to  keep  consistency  among  units. 


(4.29) 


(4.30) 
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The  evaluated  for  the  characteristics  in  Table  4.1  is 
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-0.0020 
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(4.31) 

The  eigenvalues  for  the  above  closed  loop  dynanaics  matrix  Ad  are  all  in  the  left  half  plane 
and  are  shown  in  Table  4.4.  Therefore,  the  system  is  stable  for  the  selected  gains.  It  is  also 
evident  from  the  system  responses  shown  later  in  Chapter  V  that  the  system  is  stable. 

Table  4.4  Closed  Loop  Eigenvalues 
Eigenvalues 

-0.6627  I 

-14.47 

-0.0667 

-220.37 

-1.69  +  9.06i 

-1.69  -  9.06i 

-1.29 

-0.1256 

-0.0818 


4-18 


Robustness 

A  robust  controller  is  one  that  responds  well  to  very  dynamic  plants.  It  was  stated 
in  Chapter  I  that  the  window  of  opportunity  for  a  reduction  in  drag  and  increase  in  lift 
is  small.  Since,  the  change  in  drag,  AD,  and  the  change  in  lift,  A£,  are  very  complicated 
3- dimensional  functions,  individual  slices  are  shown  for  both  the  y  and  z  directions.  The 
AD  and  AX,  are  proportional;  the  only  difference  is  the  magnitude  of  the  AX  is  about  ten 
times  greater  than  the  AD.  Therefore,  all  discussions  for  the  AD  function  also  apply  to 
the  AX  function. 

Figure  4.8  shows  a  z-slice  of  the  AD  function  at  ^  =  0.  Figure  4.9  shows  the  rate 
of  change  in  AD  evaluated  for  ^  =  0  and  y  varied  between  3.5  ft  to  43.5  ft  or  -/+§&  from 
nominal,  y  =  It  can  be  seen  from  these  two  figures  that  the  nominal  operating  point 
is  highly  dynamic.  A  slight  movement,  either  to  the  right  or  to  the  left,  can  cause  a  large 
rapid  magnitude  swing  and  even  a  sign  change  in  AD.  AD  is  almost  symmetric  with 
respect  to  y  =  for  small  changes  in  y  separation.  The  function  loses  it’s  symmetry  as 
the  deviations  from  nominal  exceed  «  —6.  This  is  easily  seen  in  Figure  4.9  where  the  left 
hand  side  forms  a  parabola  and  the  right  slowly  decays  to  zero.  Decreases  in  y  position 
greater  than  ^6  cause  the  Wing  aircraft  to  move  into  the  region  in  which  the  effect  from 
the  Lead’s  vortex  is  downwash  instead  of  upwash.  This  region  is  directly  between  the  two 
filaments  composing  the  horseshoe  vortex.  However,  the  vortex  never  becomes  a  downwash 
for  an  increase  in  y  separation  greater  than  the  nominal,  but  slowly  decays  to  zero  upwash. 

Figure  4.10  shows  the  same  AD  function,  but  sliced  in  the  y  direction  at  y  = 
Figure  4.11  shows  the  rate  of  change  in  AD  evaluated  for  y  =  and  z  varied  between 
-20  ft  to  20  ft  or  from  nominal.  It  is  obvious  from  these  two  figures  that  AD  is  also 

very  dynamic  for  changes  in  z,  but  unlike  the  z-slice  it  is  perfectly  symmetric.  The  radius 
for  a  AD  benefit  is  larger  for  the  y-slice,  about  55%  of  the  wingspan. 
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It  can  be  seen  from  the  z  and  y  slices  that  the  mciximum  reduction  in  drag  is  obtained 
at  the  nominal  y  and  z  values.  Deviations  in  the  y  and  z  positions  from  this  nominal  location 
are  compounded.  So  a  combined  change  in  y  and  z  position  in  the  formation  results  in  a 
larger  reduction  in  AD  than  for  an  individual  change  of  greater  magnitude  in  either  the  y 
or  the  z  positions. 

Figure  4.12  shows  the  actual  AD  function  and  the  first  and  second  order  models  for 
the  region  of  +/-^b  in  which  AD  is  roughly  symmetric.  It  can  be  seen,  that  the  first 
order  model  truly  does  trivialize  the  effect  of  the  upwash.  A  decrease  in  the  y  separation 
distance  causes  a  small  decrease  in  the  AD;  an  increase  in  the  y  separation  causes  a  small 
increase  in  the  AD.  However,  the  true  response  is  more  like  a  parabola.  A  decrease  in 
y  separation  causes  a  decrease  in  AD,  but  it  is  much  faster  than  the  linear  version.  An 
increase  in  the  y  separation  also  causes  a  decrease  in  the  actual  AD,  but  the  first  order 
model  causes  an  increase  in  AD  for  an  increase  in  y.  The  second  order  model  is  parabolic 
like  the  actual  and  appears  to  be  better  than  the  first  order  model  for  small  positive  and 
negative  changes  in  y,  but  for  larger  changes  the  first  order  model  is  actually  closer  to  the 
true  response  as  the  second  order  model  becomes  increasingly  negative.  The  first  order 
model  is  truly  only  accurate  for  changes  in  y  of  0.25  ft  in  magnitude.  It  can  also  be  seen 
from  Figure  4.12,  that  the  first  order  model  is  not  constant,  but  has  a  small  slope.  The 
maximum  reduction  in  drag  is  not  exactly  at  the  point  y  =  ox  23.562  ft,  but  at  23.612 
ft.  This  shift  is  caused  by  the  inclusion  of  the  correction  term  fi.  There  is  no  effect  on  the 
first  and  second  order  AD  models  for  a  change  in  z,  because  f  =  0  at  the  nominal  position. 
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Figure  4.12  A  D  at  z  =  0:  1st  Order  Approximation- Top  Actual-Middle,  2nd  Order 
Approximation-Bottom 

The  change  in  side  force,  AY  function  is  also  a  complicated  three  dimensional  surface. 
Figure  4.13  shows  a  z-slice  of  the  AY  function  at  f  =  0  and  Figure  4.14  shows  the  change 
in  side  force  for  a  change  in  y  position  evaluated  at  ^  =  0.  From  these  figures,  it  is  evident 
that  the  AY  curve  is  also  dynamic,  but  the  highly  dynamic  region  is  removed  from  the 
nominal  operating  point.  Recall,  the  nominal  separation  distances  were  choosen  to  achieve 
a  maximum  reduction  in  drag  not  side  force.  For  this  reason,  the  AY  is  less  dynamic  than 
the  AD  near  the  nominal  operating  point.  However,  the  function  is  not  symmetric  for 
equal  positive  and  negative  changes  in  y  location.  A  reduction  in  y  separation  distance 
of  more  than  results  in  a  strong  increase  in  negative  side  force.  The  side  force  slowly 
decays  to  zero  for  increasing  y  separation. 

Figure  4.15  shows  the  same  AY  function,  but  sliced  in  the  y  direction  at  y  =  ^6. 
The  partial  derivatives  for  AY  with  respect  to  a  change  in  z  are  shown  in  Figure  4.16. 
It  is  again  obvious  in  these  figures  that  the  nominal  position  is  removed  from  the  highly 
dynamic  region.  A  change  in  z  separation  will  have  only  minor  impact  on  AY.  The 
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magnitude  changes  caused  by  a  change  in  z  are  much  less  severe  than  those  caused  by  a 
change  in  y. 
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Overall,  it  was  seen  from  Figures  4.8  -  4.16  with  the  exception  of  Figure  4.12,  that  the 
vortex  strength  changes  drastically  from  location  to  location  and  changes  very  quickly,  but 
the  magnitude  of  the  cross  coupling  due  to  the  vortex  is  much  smaller  than  the  magnitude 
of  the  kinematically  induced  coupling.  This  is  easily  seen  from  the  linearized  system 
dynamics  A  matrix  in  Equation  4.18.  The  Ai  4,  A3  4,  and  A4,4  terms  are  the  kinematic 
induced  coupled  terms.  The  Ai^s,  A2,3,  A3, 3,  A4,3,  As^s,  A6,3,Ai,5,  A3, 5,  and  A4,5  terms 
are  the  aerodynamic  cross  coupling  caused  by  the  Lead  vortex.  Except  for  A6,3,  all  the 
vortex  terms  compared  to  the  kinematic  cross  coupling  terms  are  quite  small.  Therefore, 
the  kinematic  coupling  has  a  stronger  influence  on  the  robustness  of  the  controller  than 
does  the  vortex  induced  aerodynamic  cross  coupling.  The  formation  flight  controller  was 
previously  shown  to  be  robust  enough  to  handle  the  kinematic  coupling  [3].  Hence,  the 
FFCS  designed  by  Dargan  is  quite  robust  with  respect  to  the  minor  disturbances  caused 
by  the  vortex,  as  well. 

4.5  Complete  Nonlinear  Simulation  Model 

The  full  linear  system  was  developed  in  the  beginning  of  this  chapter.  This  was 
needed  for  the  design  of  the  close  formation  flight  controller.  Once  the  close  FFC  is 
developed  using  the  linear  system,  it  needs  to  be  tested  on  the  actual  nonlinear  system. 

The  nonlinear  system  is  developed  by  replacing  the  first  order  close  formation  stabil¬ 
ity  derivatives  in  Equations  (3.39)  -  (3.41)  with  the  actual  coefficients  for 
and  ACyy,  shown  in  Equations  (3.15),  (3.18),  and  (3.21)  and  adjoining  them  to  the  non¬ 
linear  kinematics  Equations  (2.22)  -  (2.24).  The  full  nonlinear  system  is: 
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where  the  vertical  separation  z  =  hw  —  hL  and  the  Wing  and  Lead  aircraft  are  assumed 
to  have  the  same  vertical  dynamics  and  Equation  (4.35)  for  ip  is  inserted  into  the  x  and 
y  separation  differential  equations  as  in  the  linear  version.  This  full  nonlinear  model  and 
the  nonlinear  model  with  the  linear  close  formation  stability  derivatives,  Equations  (4.1) 
-  (4.7)  are  used  in  the  performance  evaluations  in  Chapter  V.  The  full  nonlinear  model, 
Equations  (4.32  -  4.37),  is  shown  in  Figure  4.17. 


Vw 


Figure  4.17  Formation  Flight  Control  System  Nonlinear  Simulation  Diagram 


4.6  Summary 

In  this  chapter,  the  complete  six  state  model  was  developed  and  the  method  for  gain 
selection  was  described.  The  linear  closed  loop  system  was  developed  and  a  stability  anal¬ 
ysis  was  conducted  for  the  closed  loop  system.  The  aerodynamic  coupling  was  evaluated 
and  shown  to  be  weak  in  comparison  to  the  coupling  created  by  the  kinematics,  thereby 
justifying  the  use  the  large  formation  flight  controller.  Finally,  the  large  formation  flight 
controller  was  shown  to  be  robust  enough  to  handle  the  aerodynamic  coupling  created  by 
the  vortex. 
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V.  Performance  Evaluation 


This  chapter  describes  the  complete  nonlinear  simulation  model,  including  a  de¬ 
scription  of  the  command  prefilter.  A  description  of  the  format  for  the  response  plots  is 
provided.  Performance  evaluations  for  six  cases  are  performed:  The  cases  evaluated  are 
+/-  heading  changes,  +/-  velocity  changes,  and  +/-  altitude  changes.  Each  case  is  subdi¬ 
vided  into  three  runs.  The  system  response  for  each  case  is  discussed  and  comparisons  are 
made  between  cases. 

The  simulation  model  is  shown  in  Figure  4.17.  The  Lead  and  Wing  Aircraft  blocks 
are  represented  by  the  linear  F-16  class  aircraft  ECS  developed  in  Chapter  II.  The  FFCS 
block  is  represented  by  the  model  developed  in  Chapter  IV.  The  kinematics  block  is 
modeled  using  the  fully  non-linear  Equations  (2.22)  -  (2.24)  outlined  in  Chapter  II.  The 
Horseshoe  Vortex  block  is  actually  composed  of  two  models  described  in  Chapter  III.  One 
is  the  linear  model  developed  in  Equations  (3.35)  -  (3.38).  The  other  is  the  fully  non¬ 
linear  version  developed  in  Equations  (3.15),  (3.18),  and  (3.21).  The  coupling  terms  are 
appropriately  dimensionalized  and  then  added  as  disturbances  to  the  Wing  flight  control 
system.  Both  the  linear  cross  coupling  version  and  the  non-linear  vortex  version  are  used 
in  the  following  simulations.  A  command  prefilter  is  contained  in  the  Formation  Maneuver 
Command  block.  It  processes  the  input  step  commands  so  they  follow  more  closely  the 
commands  input  into  an  actual  aircraft  autopilot.  This  basically  ramps  the  commands  into 
the  aircraft  FCS  instead  of  inputting  an  instantaneous  step  change.  This  allows  the  errors 
to  grow  gradually  instead  of  instantly  changing  from  zero  error  to  100  percent  error  and 
prevents  the  immediate  saturation  of  the  aircraft  flight  control  system.  Saturation  of  the 
FCS  for  long  periods  can  result  in  instability.  The  complete  Matlab  Simulink  simulation 
model  is  shown  in  Appendix  A. 

The  Lead  and  Wing  aircraft  are  flying  in  a  close  diamond  formation  with  x  =  60  ft 
and  y  =  23.562/t  or  2b  and  jb,  respectively  and  z  =  0.  This  formation  was  shown 
previously  in  Figure  1.2  in  which  the  Lead  is  on  the  right  side  and  in  front  of  the  Wing. 
Typical  F-16  class  aircraft  characteristics  are  listed  in  Table  4.1  and  the  evaluated  formation 
stability  derivatives  are  listed  in  Table  4.2. 


5-1 


5.1  Case  Definition 

Each  case  consists  of  three  runs;  no  cross  coupling,  linear  cross  coupling,  and  non¬ 
linear  cross  coupling.  There  are  cases  for  both  positive  and  negative  change  of  equal 
magnitude.  Comparison  of  the  positive  and  negative  cases  show  the  difference  between 
maneuvers  in  which  the  Lead  is  moving  away  from  the  Wing  and  when  the  Lead  is  moving 
toward  the  Wing.  The  positive  and  negative  heading  plots  can  also  be  used  to  determine 
the  response  in  the  event  there  are  two  Wingmen;  one  on  each  side  of  the  Lead. 

Long  distance  flying  usually  entails  flying  from  waypoint  to  waypoint  until  the  final 
destination  is  reached.  Usually  waypoints  are  defined  by  radio  beacons  or  airports.  Be¬ 
ginning  at  the  initial  airport  of  takeoff,  the  pilot  dials  in  the  autopilots’  velocity,  heading, 
and  altitude  to  reach  the  next  way  point.  Upon  reaching  the  next  waypoint,  the  pilot 
enters  the  new  velocity,  heading,  and  altitude  values  into  the  autopilot  to  proceed  to  the 
next  waypoint  and  so  on  until  the  final  destination  is  reached.  Usually,  the  paths  followed 
are  relatively  strait.  This  means  that  large  heading  changes  between  waypoints  will  not 
be  encountered.  The  simulation  cases  were  chosen  based  on  this  scenario  and  should  be 
representative  of  typical  autopilot  changes  implemented  by  a  pilot  as  he  sets  the  autopilot 
to  proceed  to  the  next  waypoint. 

The  format  for  the  time  response  plots  and  the  Lissajous  plots  is  the  same  as  described 
in  Chapter  LV  and  is  repeated  below.  The  time  plots  contain  six  subplots.  In  the  top  three 
subplots,  the  solid  line  represents  the  Lead  aircraft  and  the  dashed  line  represents  the  Wing 
aircraft.  In  some  cases  only  a  solid  line  is  visible,  this  is  due  to  the  overlap  of  the  responses 
from  the  Lead  and  the  Wing.  In  the  bottom  three  subplots,  the  solid  line  represents  the 
nominal  separation  values  and  the  dashed  line  represents  the  actual  values.  Starting  points 
in  the  Lissajous  plots  are  represented  by  an  ‘O’  and  the  ending  points  are  represented  by 
an  ‘X’.  This  format  is  used  for  each  plot  listed  in  the  following  sections. 


5-2 


5,2  Negative  30  Degree  Heading  Change 

The  formation  is  commanded  a  heading  change  of  negative  30°  in  which  the  Lead 
aircraft  turns  into  the  Wing  aircraft.  Figures  5.1  and  5.2  show  the  time  and  Lissajous 
response  plots  for  the  linearized  system  without  cross  coupling.  Even  with  the  prefilter  in 
place  there  is  a  slight  saturation  of  the  Lead  heading  response.  This  is  not  significant  and 
more  importantly  the  Lead  is  considered  a  disturbance  so  it  truly  doesn’t  matter  whether 
it  is  saturated  or  not  for  this  analysis.  The  Wing  does  not  experience  any  saturation  in 
heading  response. 

The  Wing  heading  follows  the  Lead  heading  so  closely  that  the  two  overlap.  However, 
the  Wing’s  velocity  varies  to  follow  the  heading  change  of  the  Lead.  The  Wing  must 
decrease  it’s  velocity,  since  the  Lead  is  turning  into  it.  This  causes  the  separation  in  the 
X  direction  to  initially  decrease  as  expected.  The  y  separation  initially  increases  and  then 
undershoots  the  nominal  and  slowly  increases  and  levels  off  at  the  nominal  y  separation. 
The  formation  separation  reduces  to  a  minimum  56  ft  in  the  x  direction  or  4  ft  from 
nominal.  The  closest  the  Lead  comes  to  the  Wing  in  the  y  direction  is  20.6  ft  or  3  ft  from 
nominal.  In  fact  the  maximum  positive  and  negative  deviations  in  the  y  direction  are  3  ft, 
meeting  the  first  gain  criterion.  The  x-channel  settles  within  12  seconds  and  the  y-channel 
settles  within  60  seconds.  The  formation  spacing  is  maintained  within  approximately  3  ft 
(1/10  b)  in  both  the  x  and  y  directions  and  zero  steady  state  tracking  errors  are  achieved 
after  the  heading  change  is  complete.  The  z  separation  stays  at  zero  for  the  duration  of 
the  maneuver. 

Figures  5.3  and  5.4  are  the  response  plots  for  -30°  heading  change  with  linear  cross 
coupling.  The  most  obvious  and  most  interesting  difference  between  this  run  and  the 
previous  run  without  cross  coupling  is  the  affect  of  the  cross  coupling  on  the  z-channel. 
Specifically,  the  cross  coupling  of  the  y  separation  into  the  altitude  hold  autopilot  through 
the  Ae  3  term.  The  linear  cross  coupling  causes  an  oscillation  in  the  vertical  channel.  It  is 
a  very  small  oscillation  with  a  maximum  deviation  of  3  x  10“^  ft  and  a  minimum  deviation 
of  5  X  10“®  ft  and  is  nulled  out  by  the  PI  controller.  The  initial  increase  in  y  separation 
causes  an  increase  in  the  Wing’s  lift  causing  it  to  ascend  above  the  Lead.  It  is  apparent 
from  this  response  that  the  change  in  z  separation  follows  the  linear  first  order  AL  function 
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exactly.  The  AL  first  order  model  decreases  for  decreasing  y  separation  and  increases  for 
increasing  y  separation. 

Figures  5.5  and  5.6  are  the  response  plots  for  -30°  heading  change  with  non-linear 
cross  coupling.  The  effect  of  the  non-linear  cross  coupling  is  very  different  than  that  of 
the  linear  coupling.  The  vertical  channel  perturbations  are  less  graceful  and  greater  in 
magnitude  than  that  of  the  run  with  linear  coupling.  The  x  and  y  separation  plots  are 
also  different  than  the  previous  two  runs. 

The  minimum  x  separation  increases  from  -4  ft  to  -3.85  ft.  The  maximum  x  deviation 
increases  from  0  ft  to  0.1  ft.  Basically  the  entire  x  response  plot  just  shifted  up  by  0.1 
ft.  This  is  due  to  the  change  in  drag  from  the  change  in  y  separation.  The  linearized  AD 
function,  Figure  4.12  did  not  accurately  represent  the  true  change  in  drag  for  changes  in 
y  separations  greater  than  0.25  ft  in  magnitude  from  the  nominal  y  separation. 

The  y  separation  also  changed.  The  positive  and  negative  deviations  for  the  previous 
two  runs  were  3  ft  and  -3  ft,  respectively.  In  this  run,  the  positive  and  negative  deviations 
are  2.3  ft  and  -2.2  ft,  respectively.  This  is  a  change  of  0.7  ft  and  0.8  ft,  respectively  from 
the  previous  runs.  The  error  improved  in  both  the  positive  and  negative  directions  by 
about  the  same  amount.  Recall  from  Figure  4.13,  the  nominal  point  rests  in  the  middle  of 
the  knee  of  the  curve  where  the  slope  is  changing  slowly.  In  fact,  the  change  in  side  force 
magnitude  is  almost  equal  for  changes  in  y  of  -|-/-3  ft.  The  increase  in  y  separation  causes 
the  sideforce  to  decrease  causing  the  Wing  to  turn  away  from  the  Lead.  Since,  the  Lead 
is  turning  into  the  Wing  this  speeds  up  the  Wing’s  heading  response  leading  to  a  reduced 
error  in  the  y  separation.  The  smaller  error  means  the  controller  does  not  have  to  apply 
as  much  correction  as  in  the  previous  two  runs  resulting  in  a  smaller  undershoot  as  the 
controller  pulls  the  y  separation  toward  nominal. 

The  z  separation  decreases  initially  as  y  increases.  As  the  controller  decreases  the 
y  separation  towards  the  nominal  y,  the  z  separation  begins  to  increase.  When  the  y 
separation  crosses  the  nominal,  the  z  separation  briefly  crosses  the  z  axis  becoming  slightly 
positive  and  then  quickly  decreases  again  as  y  continues  to  decrease.  The  y  separation  then 
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slowly  increases  and  levels  off  when  it  reaches  the  nominal  y  value.  The  z  separation  slowly 
increases  overshoots  and  then  starts  to  decrease  and  level  off  at  the  nominal  z  separation. 

There  are  four  primary  pieces  to  the  z  plot.  The  first  piece  is  the  first  negative  peak. 
Since  the  Wing  is  sitting  at  the  point  of  maximum  lift  in  the  vortex,  any  deviation  causes  a 
decrease  in  AL.  This  is  reversed  from  the  lineaj  coupling  run  in  which  the  initial  increase 
in  y  separation  caused  an  increase  in  z  separation  or  an  increase  in  AL.  The  crossing  of  the 
z  axis  is  interesting.  The  response  in  this  area  demonstrates  how  quickly  the  magnitude 
of  the  AL  changes  for  relatively  small  changes  in  separation  between  aircraft.  Since  the 
changes  in  z  separation  are  so  small,  the  AD,  AL,  and  AFy  forces  are  dependent  on  the 
change  in  y  separation  only. 
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5,3  Positive  30  Degree  Heading  Change 

The  positive  30°  heading  change  without  coupling  is  shown  in  Figures  5.7  -  5.8.  This 
again  causes  a  slight  saturation  in  the  Lead  heading  response.  The  two  headings  overlap, 
however  the  velocity  increases  instead  of  decreasing  as  in  the  previous  case.  In  this  case, 
the  Lead  is  turning  away  from  the  Wing  causing  the  Wing  to  play  catch  up.  Since  the  Wing 
is  on  the  outside  it  needs  to  cover  a  longer  distance  than  the  Lead  in  order  to  maintain  the 
formation.  This  causes  the  Wing  to  briefly  increase  it’s  velocity.  The  response  in  the  x  and 
y  plots  is  exactly  the  opposite  of  the  —30°  heading  change  maneuver.  Instead  of  decreasing 
as  in  the  —30°  case,  the  x  separation  initially  increases.  The  y  separation  initially  decreases, 
then  overshoots  and  slowly  decays  to  nominal.  The  maximum  separation  in  the  x  direction 
is  63.3  ft  or  +3.3  ft  from  nominal.  It  can  be  seen  from  the  plots,  that  the  closest  the  two 
aircraft  come  to  one  another  in  the  y  direction  is  20.6  ft  and  the  farthest  is  26.5  ft,  i.e., 
approximately  -3  ft  and  +3  ft  from  nominal,  respectively.  As  expected,  there  is  no  change 
in  the  vertical  channel. 

The  positive  30°  heading  change  with  linear  coupling  is  shown  in  Figures  5.9  -  5.10. 
The  primary  difference  between  the  plots  with  linear  cross  coupling  versus  those  without 
cross  coupling  is  again  the  effect  on  the  vertical  channel.  The  cross  coupling  causes  an 
initial  oscillation  in  the  vertical  channel  that  is  almost  exactly  the  opposite  of  the  previous 
case,  -30°  heading  change  with  linear  coupling.  The  PI  controller  nulls  the  perturbation 
out  after  about  60  sec.  The  perturbation  in  the  z-channel  is  again  caused  by  the  linking 
term  Ae  3  and  responds  exactly  as  expected  from  the  AX  first  order  model. 

The  non-linear  coupling  plots  are  shown  in  Figures  5.11  -  5.12.  The  x  separation 
error  is  0.1  ft  larger  than  in  the  previous  two  runs;  3.4  ft  versus  3.3  ft.  As  in  the  -30°  case 
with  non-linear  coupling,  the  flrst  order  model  does  not  accurately  represent  the  changes 
in  magnitude  of  drag  for  changes  in  y  greater  than  +/-0.25  ft. 

In  the  case  of  the  negative  and  positive  30°  heading  changes  with  linear  cross  coupling, 
the  results  were  identical  except  flipped  in  the  z  direction.  This  is  not  true  for  the  non-linear 
cases.  The  z  separation  responses  for  both  a  positive  and  negative  30°  heading  change  are 
remarkably  similar.  The  initial  decrease  in  y  causes  a  decrease  in  z.  As  y  begins  to  increase 
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the  z  separation  begins  to  increase.  The  z  separation  briefly  crosses  its  nominal  as  y  crosses 
its  nominal  and  then  starts  to  decrease  again  as  y  continues  to  increase.  The  y  separation 
then  decreases  slowly  back  to  nominal.  The  z  separation  also  increases  back  to  nominal, 
but  overshoots  slightly.  The  first  z  crossover  is  caused  by  both  the  controller  attempting  to 
null  out  the  error  and  the  increasing  AX  caused  by  the  shrinking  error  in  the  y  separation 
between  the  Wing  and  the  Lead  aircraft.  The  reason  this  response  is  almost  identical  to 
that  of  the  negative  heading  change  is  because  the  AX  function  is  nearly  symmetric  with 
respect  to  the  nominal  y  =  23.562  ft  for  y  separations  in  the  range  of  20.5  ft  -  26.5  ft  or 
-/+  3  ft  as  previously  described  and  seen  in  Figure  4.8.  The  magnitude  of  the  negative 
deviation  is  slightly  larger  (0.06  ft)  than  in  the  -30°  run  with  nonlinear  coupling  due  to 
the  slight  non-symmetry  in  AX.  As  in  the  negative  heading  change  cases,  the  Ai),  AX, 
and  AFy  depend  on  changes  in  y  separation  only. 

Recall  from  Chapter  IV,  the  second  requirement  for  gain  selection  was  to  achieve 
consistent  results  for  both  positive  and  negative  maneuvers  of  equal  magnitude.  This  is 
best  seen  by  comparing  the  Lissajous  plots,  Figures  5.2  and  5.8  for  the  —30°  and  -f30° 
runs.  The  positive  heading  change  Lissajous  plot  is  almost  identical  as  the  the  negative 
heading  change  Lissajous  plot  except  it  is  rotated  by  180  degrees.  This  is  also  true  for 
Figures  5.4  and  5.10  and  and  Figures  5.6  and  5.12 
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5.4  50  ft/sec  Velocity  Decrease 

In  this  case,  the  formation  is  commanded  a  decrease  in  velocity  of  50  ft/sec  without 
cross  coupling  as  seen  in  5.13  -  5.14.  The  velocity  of  the  Wing  follows  the  Lead  so  closely 
that  the  two  cannot  be  distinguished  on  the  velocity  time  plot.  However,  the  x  separation 
distance  immediately  decreases.  This  shows  that  there  is  some  delay  between  the  time 
the  Lead  starts  to  decrease  it’s  velocity  and  the  time  the  Wing  starts  to  decrease  it’s 
velocity.  This  is  expected  since  the  Wing  does  not  have  foreknowledge  of  the  Lead’s  future 
maneuvers.  The  FFCS  is  very  tight  and  allows  a  maximum  error  of  approximately  0.8 
ft  and  settles  in  about  60  seconds.  The  y  and  z  channels  are  unaffected  by  the  velocity 
change. 

There  is  no  difference  between  the  run  with  cross  coupling  and  the  run  with  no  cross 
coupling.  As  seen  in  Figures  5.15  -  5.16,  the  response  in  the  x  separation  is  the  same 
as  the  previous  run  without  cross  coupling.  The  first  order  model  was  developed  based 
on  the  assumption  that  the  vortex  extends  a  long  distance  behind  the  Lead.  Secondly, 
the  linearized  kinematic  equation  for  y  separation  is  only  affected  by  heading  changes,  not 
velocity  changes.  Therefore,  a  change  in  x  direction  does  not  affect  the  upwash  or  sidewash. 
Only  a  change  in  the  y  direction  can  cause  a  change  in  upwash  and  sidewash.  Therefore, 
the  reduction  in  drag  AD  is  constant  for  pure  velocity  changes.  Prior  to  the  Lead  velocity 
change  in  this  run  and  the  previous  run  the  Wing  aircraft  was  in  steady-state,  i.e.,  aU  the 
forces  in  the  x  direction  sum  to  zero,  leading  to  a  constant  velocity.  Therefore,  the  error 
caused  by  the  decrease  in  the  Lead’s  velocity  and  input  into  the  x-channel  of  the  FFC  is 
of  equal  magnitude  and  sign  and  causes  identical  responses  in  both  runs. 

The  responses  to  a  50  ft/sec  decrease  in  velocity  with  non-linear  cross  coupling  can 
be  seen  in  Figures  5.17  -  5.18.  This  is  another  very  interesting  run.  The  x  response  is 
the  same  as  in  the  previous  two  runs.  However,  the  heading,  altitude,  y  separation,  and  z 
separation  responses  are  not  zero  as  in  the  previous  two  runs.  Both  the  Wing’s  heading 
and  y  separation  plots  experience  an  initial  oscillation.  The  vertical  channel  experiences  a 
slight  decrease  in  AX  and  causes  a  drop  in  altitude  of  2  x  10“®  ft.  The  inital  decrease  in  y 
separation  causes  a  slight  decrease  in  AX  that  leads  to  an  initial  decrease  in  z  separation 
and  then  a  slow  return  to  nominal.  Equations  (3.15),  (3.18),  and  (3.21)  are  necessary  to 
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understand  the  responses  from  this  run.  The  linear  model  assumes  that  the  formation’s 
velocity  is  constant  at  825  ft/sec.  This  is  not  significant  for  the  heading  change  maneuvers, 
since  the  Lead  and  Wing  velocities  vary  only  slightly.  However,  this  is  not  the  case  for 
a  velocity  change  of  50  ft /sec.  There  are  two  reasons  for  this  outcome.  First,  the  Lead 
velocity  determines  the  strength  of  the  vortex.  Secondly,  the  Wing’s  velocity  determines 
the  change  in  angle  of  attack,  a,  due  to  the  upwash.  Both  of  these  directly  affect  the 
strength  of  the  upwash  and  downwash  exerted  on  the  Wing  aircraft.  The  AX  function 
is  proportional  to  the  Lead’s  velocity  and  inversly  proportional  to  the  Wing’s  velocity,  as 
seen  in  Equation  (3.18).  The  formation  hold  controller  responds  very  quickly  keeping  the 
velocity  of  the  Wing  almost  equal  to  the  velocity  of  the  Lead.  Since  the  two  velocities  are 
almost  equal,  the  change  in  angle  of  attack,  Aa,  is  almost  constant.  It  decreases  the  most 
at  the  beginning  of  the  maneuver  when  the  error  between  the  Lead  and  Wing  velocities  is 
the  greatest.  This  is  seen  in  the  initial  decrease  in  z  separation.  AXy  is  affected  in  exactly 
the  same  way  as  AD  and  AX.  This  leads  to  the  tear  drop  shape  shown  in  Figure  5.18. 
Even  though  this  response  looks  very  dynamic,  closer  inspection  reveals  the  magnitudes 
of  the  y  deviations  are  very  small  and  the  x  separation  is  essentially  unchanged  from  the 
previous  two  runs. 
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50  ft /sec  Velocity  Decrease  With  Linear  Coupling  Lissajous  Plot 
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Figure  5.17  50  ft/sec  Velocity  Decrease  With  Non-Linear  Coupling  Time  Plot 


Figure  5.18  50  ft/sec  Velocity  Decrease  With  Non-Linear  Coupling  Lissajous  Plot 


5,5  50  ft/ sec  Velocity  Increase 

The  formation  is  commanded  an  increase  in  velocity  of  50  ft/sec  without  cross  cou¬ 
pling.  This  case  is  the  exact  opposite  of  the  previous  case  run,  50  ft/sec  decrease  in  velocity 
with  no  coupling.  The  velocity  of  the  Wing  again  follows  the  Lead  so  closely  that  the  two 
cannot  be  distinguished  on  the  time  plot,  Figure  5.19,  but  the  x  separation  distance  imme¬ 
diately  increases  as  seen  in  the  Lissajous  plot.  Figure  5.20.  As  in  the  previous  case,  there 
is  some  delay  between  the  time  the  Lead  starts  to  increase  velocity  and  the  time  the  Wing 
starts  to  increase  velocity.  The  FFCS  allows  a  maximum  error  of  1.4  ft,  0.6  ft  greater  than 
the  previous  case  and  settles  in  about  60  seconds.  Figures  5.19  and  5.20.  The  difference  in 
magnitudes  between  cases  is  caused  by  the  rate  limits  for  the  Mach  hold  autopilot.  The 
upper  limit  is  half  of  the  lower  limit.  Therefore,  the  aircraft  can  decelerate  faster  than 
they  can  accelerate.  The  point  at  which  the  Wing’s  Mach  hold  autopilot  saturates  is  easily 
identified  by  the  constant  slope  of  the  x  separation  plot  for  the  first  6  seconds.  Except 
for  the  affects  of  the  non- symmetric  saturation  limits,  the  second  gain  criterion  is  again 
satisfied  for  this  case  and  the  case  of  50  ft /sec  velocity  decrease. 

Figures  5.21  -  5.22  show  the  aircraft  response  to  linear  cross  coupling.  This  is  identical 
to  the  previous  run  without  coupling.  Recall  from  the  previous  case  for  a  velocity  decrease 
of  50  ft /sec;  the  linear  first  order  model  for  lift,  drag,  and  sideforce  is  not  affected  by 
changes  in  x  separation  and  velocity. 

The  responses  for  an  increase  in  velocity  of  50  ft /  sec  with  non-linear  cross  coupling 
is  shown  in  Figures  5.23  -  5.24.  As  in  the  case  with  non-linear  cross  coupling  for  the 
velocity  decrease  of  50  ft /sec  case,  all  channels  are  affected  by  the  velocity  change.  The  y 
separation  experiences  an  initial  small  oscillation.  The  x  separation  error  is  the  same  as 
the  two  50  ft/sec  velocity  increase  runs  without  coupling  and  with  linear  coupling.  The 
vertical  channel  experiences  an  initial  increase  in  lift  caused  by  the  initial  increase  in  Aa 
caused  by  the  Lead’s  velocity  increasing  faster  than  the  Wing’s  velocity.  Rotating  this 
Lissajous  plot  and  comparing  with  the  Lissajous  Plot  from  the  -50  ft/sec  run.  Figure  5.18 
reveals  the  responses  are  very  similar  except  in  this  run,  the  y  deviation  is  shifted  left  of 
the  nominal.  This  results  in  the  nonsymmetrical  bulge  in  the  inverted  tear  drop. 
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Time  (sec) 


Figure  5.24  50  ft/sec  Velocity  Increase  With  Non-Linear  Coupling  Lissajous  Plot 


5.6  1000  ft  Altitude  Decrease 

In  this  case,  the  formation  is  commanded  a  1000  ft  decrease  in  altitude.  Figures  5.25  - 
5.26  show  the  response  without  cross  coupling.  The  altitude  channel  is  completely  separate 
from  the  horizontal  channels  when  cross  coupling  is  not  included.  The  Wing  tracks  fairly 
closely,  but  does  deviate  slightly  from  the  Lead.  This  can  be  seen  in  the  z  time  plot, 
which  shows  the  vertical  separation  grow  very  quickly,  level  off  for  a  few  seconds  and  then 
begin  to  decrease.  This  maneuver  saturates  both  the  Lead  and  Wing  altitude  autopilots. 
The  PI  controller  experiences  windup  during  this  period  of  saturation.  This  can  be  seen 
from  the  level  area  between  1  and  5  seconds  in  the  z  separation  plot.  Integrator  windup 
causes  the  feedback  to  be  much  larger  than  the  command  leading  to  a  large  error  fed  into 
the  controller  even  though  the  actual  response  may  be  close  to  the  commanded  response. 
In  this  case,  the  windup  only  causes  a  delay,  but  in  extreme  circumstances  windup  can 
cause  the  aircraft  to  become  unstable.  Prefilters  are  used  to  prevent  this  from  occurring 
in  most  cases.  They  slowly  ramp  the  command  into  the  system,  allowing  the  error  to 
grow  gradually  and  give  the  controller  time  to  respond.  The  prefilter,  described  at  the 
beginning  of  this  chapter,  is  working,  but  not  as  well  as  desired.  The  aircraft  does  not  go 
unstable,  but  does  take  a  few  seconds  to  winddown.  Figure  5.26  shows  there  is  no  affect 
on  the  horizontal  separation  distances  from  the  change  in  altitude. 

The  responses  for  a  1000  ft  decrease  in  altitude  with  linear  cross  coupling  are  shown 
in  Figures  5.27  -  5.28.  Recall  from  Sections  5.3  and  5.4,  a  change  in  heading  with  linear 
cross  coupling  caused  a  perturbation  in  the  vertical  channel.  This  case  shows  the  reverse  is 
true  also.  A  change  in  altitude,  with  cross  coupling  included,  leads  to  a  minor  oscillation 
in  the  heading.  The  A4^5  term  directly  links  the  z  separation  into  the  heading  autopilot. 
This  is  easily  seen  in  the  plot  for  the  heading  response.  The  heading  is  linked  into  the 
Mach  hold  autopilot  and  the  x  and  y  separations.  This  causes  the  small  initial  oscillations 
in  the  velocity  and  the  x  and  y  separation  plots.  The  maximum  y  separation  magnitude 
is  2  X  10~®  ft.  Under  real  circumstances  all  of  these  effects  would  be  drowned  out  by 
insignificant  changes  in  air  currents,  but  it  is  still  interesting  to  see  the  simulation  behave 
as  expected.  Overall,  everything  settles  down  in  about  50  seconds  with  the  exception  of 
a  little  bit  of  chattering.  The  linear  first  order  model  AL  and  AD  is  not  affected  by  a 
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change  in  z.  All  the  changes  in  AL  and  AD  are  caused  by  the  change  in  heading  which  is 
caused  by  a  change  in  sideforce,  AFy^  due  primarily  to  the  change  in  z  separation,  since 
the  y  separation  is  extremely  small.  Therefore,  the  affect  of  a  change  in  z  separation  on 
the  X“Channel  is  only  superficial  and  the  affect  on  the  y-channel  is  minor  due  to  the  simple 
first  order  model  employed. 

Figures  5.29  -  5.30  show  the  response  for  a  decrease  of  1000  ft  in  altitude  with  non¬ 
linear  cross  coupling.  The  inter-relationship  between  the  z  and  y  separations  is  much  more 
obvious  for  this  run  than  for  the  linearized  run.  The  comparison  between  this  run  and  the 
previous  shows  the  x  and  y  separation  errors  are  much  larger,  highlighting  the  difference 
between  the  full  non-linear  model  and  the  model  with  linear  aerodynamic  cross  coupling. 
The  y  separation  has  a  maximum  error  of  0.5  ft  and  a  minimum  error  of  -0.2  ft.  The 
maximum  x  separation  error  is  also  much  larger  in  magnitude,  0.1  ft.  The  y  deviations 
are  relatively  small  compared  to  the  heading  change  cases,  but  larger  than  those  from  the 
velocity  maneuver  change  cases.  Recall  from  Figure  4.13,  the  change  in  sideforce,  AFy, 
was  not  as  dynamic  as  the  change  in  drag  for  small  changes  in  y  separation.  However,  as 
seen  in  Figure  4.15  large  changes  in  z  separation  can  cause  large  changes  in  sideforce.  A  z 
separation  of  6  ft  leads  to  an  increase  in  magnitude  of  AFy  of  50%.  Since  AFy  is  negative 
it  pushes  the  Wing  away  from  the  Lead,  increasing  the  y  separation.  The  y  separation 
has  a  maximum  magnitude  of  1.1  ft.  Since  the  nominal  point  sits  on  the  knee  of  the  curve 
as  seen  in  Figure  4.13,  AFy  does  not  vary  much  for  the  1.1  ft  change  in  y  separation. 
Therefore,  in  this  case  the  effects  of  the  change  in  z  separation  are  greater  than  the  effects 
from  the  change  in  y  separation. 
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5»  7  1000  ft  Altitude  Increase 

The  formation  is  commanded  a  1000  ft  increase  in  altitude  without  coupling  as  shown 
in  Figures  5.31  -  5.32.  This  case  is  almost  identical  to  the  previous  except  the  direction 
for  each  term  is  opposite.  The  Wing  tracks  fairly  closely,  but  does  deviate  from  the  Lead. 
This  can  be  seen  in  the  z  time  plot,  which  shows  the  vertical  separation  grows  very  quickly, 
levels  off  for  a  few  seconds  and  then  decreases.  This  maneuver  saturates  both  the  Lead  and 
Wing  altitude  autopilots,  but  the  Wing  is  saturated  longer  than  the  Lead.  This  is  made 
better  and  worse  by  the  upper  limit  on  the  FCS  altitude  hold  autopilot.  The  magnitude 
of  the  upper  limit  is  1/3  less  than  the  lower  limit.  The  smaller  limit  causes  the  Lead 
autopilot  to  saturate  faster,  thereby  not  allowing  the  Lead  to  move  farther  away  from  the 
Wing.  However,  the  lower  saturation  limit  causes  the  vertical  PI  controller  to  windup  for 
a  longer  period  of  time,  causing  a  longer  delay  before  the  controller  winds  down  and  starts 
bringing  the  z  separation  back  to  the  nominal  value.  This  is  easily  seen  in  the  z  separation 
in  which  it  is  constant  from  1  to  7  seconds,  two  seconds  longer  than  in  the  1000  ft  altitude 
decrease  case.  There  is  no  effect  on  the  horizontal  separation  distances  from  the  change 
in  altitude.  The  reasons  for  this  were  described  in  the  last  section  in  the  case  of  a  1000  ft 
altitude  decrease. 

Linear  cross  coupling  causes  some  interesting  results  as  in  the  previous  case.  The 
responses  for  this  run  are  shown  in  Figures  5.33  and  5.34.  A  change  in  altitude  with 
cross  coupling  included  initially  causes  minor  oscillations  in  the  heading  and  velocity.  The 
changes  in  heading  and  velocity  lead  to  slight  perturbations  in  the  x  and  y  separation 
distances.  Overall,  everything  settles  down  in  about  60  seconds.  The  x  separation  continues 
to  oscillate,  but  the  magnitude  of  the  oscillation  is  on  the  order  of  nanofeet.  The  effects 
on  the  horizontal  channels  is  again  due  to  the  coupling  caused  by  the  term. 

Figures  5.35  -  5.36  show  the  response  for  an  increase  of  1000  ft  in  altitude  with  non¬ 
linear  cross  coupling.  As  in  the  previous  -1000  ft  case  with  nonlinear  cross  coupling,  the 
velocity  and  heading  responses  have  an  initial  small  oscillation.  The  y  separation  has  a 
minimum  error  of  -1.8  ft  and  a  maximum  error  of  0.5  ft.  Recall  from  the  previous  case, 
a  z  error  of  6  ft  resulted  in  an  increase  in  sideforce  of  50%.  In  this  case,  a  z  error  of  -6 
ft  results  in  a  decrease  of  130%  causing  the  AFy  to  actually  become  slightly  positive  as 
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seen  in  Figure  4.15.  This  strong  positive  AFy  pulls  the  Wing  toward  the  Lead  causing 
the  y  separation  to  decrease.  Since  the  AFy  is  stronger  for  a  negative  z  change  than  for  a 
positive  z  change  this  causes  the  maximum  magnitude  in  y  separation  to  be  0.7  ft  greater 
between  this  run  and  the  -1000  ft  altitude  with  nonlinear  coupling  run.  As  expected,  the 
maximum  x  separation  error  is  larger  in  magnitude,  0.1  ft,  than  the  linear  run.  Overall, 
all  perturbations  settle  in  about  60  seconds. 

It  is  again  apparent  from  Figures  5.25  -  5.36  that  the  results  are  consistent.  The 
error  in  the  z  channel  is  about  6  ft  for  both  an  increase  and  decrease  in  altitude.  The  x 
and  y  changes  also  show  some  consistency  between  responses  for  an  increase  and  decrease 
in  altitude.  Therefore,  these  results  satisfy  the  second  gain  criteria  of  consistent  results  for 
positive  and  negative  maneuver  changes. 
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5. 8  Summary 

Six  cases  were  presented:  +/-30®  heading  change,  +/-50  ft/sec  velocity  change,  and 
+/-1000  ft  altitude  change.  Some  observations  not  mentioned  in  the  cases  are  presented: 
The  FFC  can  maintain  the  formation  within  the  10%  tolerance  for  y  and  z  errors  for 
a  Lead  heading  change  of  +/-30°,  +/-50  ft/sec  velocity  change,  and  +/-400  ft  altitude 
change.  The  400  ft  altitude  response  was  not  shown  in  the  cases.  The  +/-30°  did  not 
meet  the  10%b  tolerance  for  x  error,  but  since  the  change  in  x  is  negligible  this  is  not  a 
requirement.  However,  a  +/-20°  heading  change  does  meet  the  10%b  for  x,  y,  and  z  errors. 
The  responses  were  also  shown  to  meet  the  gain  criteria  listed  in  Chapter  IV. 
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VL  Conclusion 


6,1  Conclusions 

The  primary  objective  of  this  research  was  to  determine  if  a  large  formation  flight 
controller  could  be  redesigned  and  applied  to  a  close  formation.  The  close  formation  flight 
controller  designed  in  Chapter  IV  has  met  this  objective. 

The  primary  question  to  be  answered  by  this  research  was  to  determine  if  a  close 
formation  flight  controller  designed  without  due  consideration  of  the  close  formation  aero¬ 
dynamic  coupling  can  handle  the  additional  forces  created  by  the  close  formation.  The 
secondary  question  was  to  determine  if  the  kinematic  coupling  is  stronger  than  the  aero¬ 
dynamic  cross  coupling.  The  primary  and  secondary  questions  have  also  been  answered 
through  the  stability  and  robustness  analysis  of  Chapter  IV  and  the  results  obtained  in 
Chapter  V. 

The  effects  of  close  formation  flight  were  derived  and  applied  to  the  close  formation 
flight  controller.  It  was  shown  that  the  large  formation  flight  controller,  with  the  gains 
adjusted  for  a  close  formation,  resulted  in  zero  steady  state  tracking  errors.  The  maximum 
tracking  error  in  x  separation  was  4  ft  («  ^6),  the  maximum  tracking  error  in  the  y 
separation  was  3  ft  and  the  maximum  tracking  error  in  the  z  separation  was  6  ft 

(|6)  for  the  envelope  defined  in  the  assumptions  with  the  full  non-linear  kinematics  and 
aerodynamic  cross  coupling  included. 

The  criteria  for  gain  selection  was  also  met.  The  gain  criteria  are  briefly  summarized 
as:  similar  positive  and  negative  perturbations  for  a  heading  change  maneuver,  similar 
results  for  a  positive  and  negative  maneuver  change,  and  the  Wing  should  not  cross  over 
the  X  or  y  axes.  Consistent  positive  and  negative  deviations  from  nominal  for  a  heading 
change  maneuver  satisfied  the  first  gain  criterion.  The  responses  for  a  negative  and  positive 
heading  and  velocity  maneuver  change  of  equal  magnitude  were  consistent;  satisfying  the 
second  gain  criterion.  Finally,  the  Wing  aircraft  never  crossed  over  either  the  x  or  y  axes. 
Therefore,  no  collisions  occurred  and  the  third  gain  criterion  was  satisfied. 

With  the  exception  of  the  Aq^s  term,  the  kinematic  coupling  was  generally  at  least 
two  orders  of  magnitude  greater  than  any  aerodynamic  induced  cross  coupling.  The  close 
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formation  flight  controller  developed  without  including  aerodynamic  cross  coupling  proved 
to  be  stable  and  robust.  All  eigenvalues  for  the  closed  loop  system  were  shown  to  be  in 
the  left  half  plane.  The  overall  responses  from  the  system  with  the  inclusion  of  the  full 
nonlinear  kinematics  and  aerodynamic  effects  were  shown  to  be  similar  to  the  responses 
from  the  full  nonlinear  kinematics  without  the  inclusion  of  the  aerodynamic  effects. 

Finally,  it  was  shown  that  the  formation  flight  controller  can  maintain  the  formation 
geometry  within  the  required  10%  tolerance  for  y  and  z  errors  for  a  Lead  heading  change 
of  +/-30°,  a  velocity  change  of  +/-50  ft/sec,  and  an  altitude  change  of  +/-400  ft.  Thus 
the  close  formation  flight  controller  can  enable  aircraft  to  take  advantage  of  the  reduction 
in  induced  drag  brought  about  by  the  aerodynamic  coupling  effects. 

6.2  Recommendations  for  Further  Study 

•  The  models  employed  in  this  research  were  rather  rudimentary.  The  prefllter  was 
a  simple  first  order  lag  and  was  identical  for  each  channel.  The  FCS  models  were 
based  on  a  point  mass  and  did  not  include  moment  effects.  The  FCS  C-130  models 
were  modified  to  represent  an  F-16  class  aircraft.  More  accurate  Models  should  be 
developed  to  better  represent  a  true  F-16  class  aircraft. 

•  Moment  effects  caused  by  the  vortex  should  be  developed.  The  effects  of  linearized 
and  non-linear  moments  should  be  applied  to  the  close  FFC  to  determine  if  it  is 
robust  enough  or  whether  it  needs  to  be  modified. 

•  If  higher  precision  formation  flight  control  is  desired,  it  may  be  prudent  to  add  sensor 
noise  to  the  system  and  perform  a  stochastic  based  analysis  using  a  Kalman  filter 
to  estimate  the  true  location  of  the  Lead  with  respect  to  the  Wing  in  the  formation. 
This  will  also  promote  fusing  sensor  data  from  a  close  formation  of  UAVs  performing 
a  surveillance  mission  in  which  each  UAV  has  a  different  sensors. 

•  Implement  the  close  formation  flight  controller  in  a  digital  controller. 

•  Incorporate  anti-collision  circuitry  and  analyze  the  effect. 
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6.3  Summary 

The  close  formation  aerodynamic  coupling  effects  on  the  Wing  aircraft  caused  by  the 
Lead’s  wing  vortex  have  been  included  in  the  formation  flight  control  system  dynamics. 
A  close  formation  flight  control  system  for  the  Wing  aircraft  was  designed.  It  was  shown 
that  formation  flight  control  system  designs,  accounting  for  kinematic  coupling  effects 
only,  are  robust  and  can  handle  the  additional  aerodynamic  coupling  effects  caused  by  close 
formation  flying.  The  developed  close  formation  flight  controller  can  enable  aircraft  to  take 
advantage  of  benefits  from  the  vortices  created  by  the  Lead  aircraft  in  the  formation.  This 
can  reduce  the  formations  fuel  consumption  and  extend  formation  range  and  endurance. 
This  can  also  enable  multiple  aircraft  or  UAVs  to  fly  together  in  a  controlled  close  formation 
during 
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Appendix  A.  Appendix  A 


The  F-16  characteristic  data  are  contained  in  the  Matlab  data  file,  controlfile.m.  The 
actual  program  is  contained  in  the  Matlab  script,  runsimc.m.  The  script  loads  the  data  file 
and  then  executes  the  simulations  using  the  Simulink  models.  Two  models  are  used  the 
first  is  the  linear  state  space  model  and  the  second  is  the  complete  linear /nonlinear  block 
model.  The  linear  simulation  is  used  to  perform  initial  analysis  and  design  of  the  controller. 

This  includes  the  gain  selection  described  in  Chapter  IV.  The  full  block  simulation  model 
is  used  to  test  the  designed  controller  for  stability  and  robustness  and  determine  if  the 
gains  need  to  be  adjusted. 

A.l  Characteristic  Data  file 

*/,  controlfile.m 
•/.  PLEASE  NOTE 

*/,  Second  order  MODELS  with  saturations  most  accurately  follow  true  autopilots 
*/,  without  saturations  if  the  inputs  are  limited  to  30  degree  magnitude  heading 
'/,  change;  40  ft/sec  magnitude  velocity  change;  1000  ft  magnitude  altitude  change 
*/,  More  thorough  methods  will  have  to  be  employed  if  values  greater  than  these 
*/,  are  required.  They  might  include  ramping  of  input  or  alteration  of  models 
'/,  based  on  better  modeling  methods  as  outlined  in  Russ  Miller’s  dissertation 

'/.  NON-LINEAR  or  LINEAR 

linear=l ; 
nonlin=-l ; 

*/, switch  =  linear 

•/,  THIS  IS  FIRST  ORDER  MODELS 
•/.  COMMANDS 

y,  F-16  wingspan  30  ft 
•/,  Pi/4  X  30  =  23.6  ft 
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nominal. vel=825 ;  '/,  M=0.85  a=971 
command_vel=825 ; 
nominal_heading=0 ; 
commanded.headlng-O ; 
nominal.alt  =  45000; 
commanded.alt  =  46000; 

x_nominal=  60;  */,  2b  nominal  this  one  does  not  matter  as  much  for  upwash 
x_ command*  60; 

y_nominal=  23.562;  */,  (Pi/4  X  132) /2  nominal  this  one,  we  also  need  a 

y. command*  23.562; 

z. nominal  =  0; 
z.command  =  0; 

y,  F-16  LEAD  AIRCRAFT  AUTOPILOT 
'/,  time  constants 
lead2_vel_time=5 ; 
leadl_Tsi_time=l/3 ; 
lead2_alta_time=0 . 3075 ; 
lead2_altb_time=3 . 85 ; 

y,  Lead  Saturations 
lead_accel_sat_low=-10; 
lead_accel_sat_high=5 ; 
lead_Tsi_sat_low=-6; 
lead_Tsi_sat_high*6; 
lead_alt_sat_low=-126 ; 
lead_alt_sat_high=100 ; 

y.  F-16  WING  AIRCRAFT  AUTOPILOT 


as  does 


z  component 
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*/,  Time  constants 
wing2_vel_timG=5 ; 

wingl_Tsi_time=l/3;  '/.for  a  true  second  order  model 
wing2_alta_time=0 . 3075 ; 
wing2_altb_t ime=3 . 85 ; 

%  Wing  Saturations 
wing_accel_sat_low=-10 ; 
wing_accel_sat_high=5 ; 
wing_Tsi_sat_low=-6 ; 
w ing_Ts i _  s  at  _high=6 ; 
wing_alt_sat_low=-126 ; 
wing_alt_sat_high=100 ; 

'/,  Stability  Derivatives  and  aircraft  values  F-16  class 
weight _max=25000 ; 

weight_avg=21500;  */,  lbs  half  tank  of  fuel 

we  ight_min=  18000;  */,  no  fuel 

AR=3;  */,  aspect  ratio 

S=300;  */,  sq  ft  area  of  wing 

b=30;  */,  ft  wingspan 

Svt=54.75;  */,  sq  ft 

hz=10  ;  */,  ft  or  120  inches  height  of  tail 
aw=5.3  ;  ’/,  per  rad  wing  lift  curve  slope 
avt=5.3  ;  '/,  per  rad  tail  lift  curve  slope 
eta=.95  ;  '/,  efficiency  factor  of  tail 
q=155.8;  '/,  Ib/sq  ft  dynamic  pressure 
mu=0.03;  ’/,  upwash/sidewash  correction  factor 
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•/.  DARGAN  CONTROLLER 


•/.  LINEAR  MIXER  cur  DAR 
vel_error_gain=12.5;  */,  5 
hGad_error_gain=6j  */,  10 
x_pos_error_gain=-8;  '/,  2 
y_pos_error_gain=-0.6;  '/,  1 
z_pos_error_gaiii=25;  */,  0.5 

'/,  PI  CONTROLLER 
x_pro_gain=6;  */.  0.17  1.2 
x_int.gain=0.4;  */,  0.02 
y-Pro_gain=ll;  */,  0.5 
y_int_gain=0.9;  */.  0.05 
2_pro_gain=4.0;  */,  1.0 
2_int_gaiii=0.5;  */.  0.5 
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A, 2  Executable  Program 

clear 

controlf ile 

def _time=60; 
def_cont=l; 
f  ig_niim=0; 
def.par=[]  ; 

y,  close  old  simulink  controller  if  open 

parms=get.param; 

if  parms  ==  def^par; 

9 

else 

close^system; 

end 

•/simplify  names  of  time  constants 
tv=lead2_vel_time ; 
tpsi=leadl_Tsi_time ; 
tha=lead2.alta.time ; 
thb=lead2_altb_time ; 

•/  simplify  names  of  gains 
Kv=vel_error_gain; 

Kpsi=head.error.gain; 
Kx=x_pos_error_gain; 
Ky=y_pos_error_gain ; 

Kz=z_pos .error. gain ; 
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y.  PI  CONTROLLER 
KXP=x_pro_gain ; 
KXI=x_int_gain; 
KYP=y_pro_gain ; 
KYI=y_int_gain; 
KZP=z_pro_gain ; 
KZI=z_int_gain; 


prev=input ( ’ Do  you  want  to  close  previous  pictures  1  yes  0  no:  ’) 

if  prev  ==  1 

close  all; 

hx=l; 

hy=0; 

else  prev  ==  0 
prev=0 ; 
hx=l ; 
liy=l; 

end 

cont=input( ’please  choose  1  for  current  or  0  for  previous  :  ’); 

if  cont  ==  1 

cd  . . / dargan/current ; 

XI=[0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0]  ; 
else  cont  ==  0 
cd  . . /dargein/old ; 

y,XI=[0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0]  ; 

XI= [0,0, 0,0, 60, 23. 562, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0] ; 
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choice-input ( ’pleas6  choose  1  for  linear  or  0  for  nonlinear  kinematics:  ’); 

if  choice  ®=  1; 

kswitch  =  linear; 

else  choice  ==  0; 

kswitch  =  nonlin; 

end 


choice=input( ’please  choose  1  for  linear  or  0  for  nonlinear  disturbance:  ’); 

if  choice  ==  1; 

dswitch  =  linear; 

else  choice  ==  0; 

dswitch  =  nonlin; 

end 


V_dist_sw=input( ’Do  you  want  a  velocity  disturbance,  1  yes  0  no:  ’); 
H_dist_sw=input(’Do  you  want  a  altitude  disturbance,  1  yes  0  no:  ’); 
Tsi_dist_sw=input(’Do  you  want  a  heading  disturbance,  1  yes  0  no:  ’); 


end 


’Calculate  formation  stability  derivatives’ 
weight=weight_max ; 

Cl=weight/(q*S) ; 
m=weight/32.2  ;  '/.mass 
bp*(pi*b)/4; 

•/.  TRUE  CHANGE  IN  DRAG 
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'/,  nominal 

G=(2*Cl*nominal_vel*b)/(pi*AR) ; 

W=G/(4*pi*bp)*(log((y_nominal“2  +  z_nominal~2  +  (b*mu)~2)/((y_nominal-bp)~2  +... 
z_nominal~2  +  (b*mu)~2)) . . . 

-  log(((y_nominal+bp)~2  +  z_nominal“2  +  (b*mu)"2)/(y_nominal'‘2  +... 
z_nominal~2  +  (b*mu)'2))); 
delta_C_D=Cl*W/nominal_vel ; 

D_nominal=q*S*delta_C_D/m 


•/,  TRUE  CHANGE  IN  LIFT 
y.  nominal 

delta_C_D=aw*W/nominal_vel; 

L_nominal=q*S*delta_C_D/m 


'/,  TRUE  CHANGE  IN  SIDE  FORCE 
'/.nominal 

V=G/(4*pi*hz)*(log(((y_nominal-bp/2)“2  +  z_nominal“2  +  (b*mu)“2)/((y_nominal-bp/2)“2  +. 
(z_nominal+hz)“2  +  (b*mu)“2)) . . . 

-  log(((y_nominal+bp/2)'‘2  +  z_nominal'‘2  +  (b*mu)“2)/((y_nominal+bp/2)“2  +... 
(z_nominal+hz)'‘2  +  (b*mu)“2))); 
delta_C_Y=eta*Svt*avt*V/(S*nominal_vel) ; 

Tsi_nominal=q*S*delta_C_Y/m 

*/,  velocity 

delta_C_D_W_y= (1/ (pi*AR) ) * (Cl~2) * ( 1/b) * ( (3/ 8) *pi/ [( (pi/4) ~2+mu''2) * ( (pi/2) “2+mu“2) ] ) ; 

'/•  altitude 

delta_C_L_W_y=aw*Cl* (3/8) *pi/  [pi>KAR*b* ( (pi/4) “2+mu~2) * ( (pi/2) “2+mu“2)] ; 
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'/•  heading 

C_Y_W_y_f irst=hz~2/ [2*b“2*((pi/8)“2+mu“2)*( (pi/8)''2+(hz/b)~2+mu~2)] ; 

C_Y_W_y_sec=3*hz“2/ [2*b“2* ( (3*pi/8) “2+mu“2) * ( (3*pi/8) “2+ (hz/b) "2+mu“2)] ; 
C_Y_W_y_inner=C_Y_W_y_f irst-C_Y_W_y_sec; 
delta_C_Y_W_y=Gta*Svt*avt*Cl*C_Y_W_y_inner/ (pi*AR*2*S*hz) ; 

delta_C_Y_W_z=- (eta*Svt*avt*Cl*32*hz*pi“2) / (pi*AR*2*S*pi*b~2*8“2* . . . 

C(  (pi/8)  “2+mu“2+  (hz/b)  ~2)  *  (  (3*pi/8)  “2+mu''2+  (hz/b)  “2) ]  )  ; 

if  cont  ==  1 

lindis=input( ’please  choose  1  for  disturb  or  0  for  nodisturb  :  ’); 

if  lindis  ==  1 

'/.  STATE  SPACE  MODELS 

*/,  A  with  disturbances 

A=[0  -1  [(q*S/(m*nominal_vel))*(delta_C_Y_W_y*y_nominal)]  [-(y_nominal/wingl_Tsi_time)*(p 
0  [-(l/wing2_vel_time)]  [(q*S/m)*(delta_C_D_W_y)]  000; 

0  0  [-(q*S/(m*nominal_vel))*(delta_C_Y_W_y*x_nominal)]  [((x_nominal/wingl_Tsi_time  -  nomi 
0  0  C(q*S/(m*nominal_vel))*(delta_C_Y_W_y)*(180/pi)]  [-(l/wingl_Tsi_time)]  [(q*S/(m*nomin 
000001; 

0  0  C(q*S/m)*(delta_C_L_W_y)]  0  [-l/(wing2_alta_time*wing2_altb_time)]  - [(l/wing2_alta_ti 

else  lindis  ==  0 

y,  A  without  disturbances 

A=[0  -1  0  [-(y_nominal/wingl_Tsi_time)*(pi/l80)]  0  0; 

0  C-(l/wing2_vel_time)]  0000; 

000  C((x_nominal/wingl_Tsi_time  -  nominal_vel))*(pi/180)]  0  0; 

000  [-(l/wingl_Tsi_time)]  0  0; 

000001; 

0000  [-l/(wing2_alta_time*wing2_altb_time)]  -[(l/wing2_alta_time)  +  (l/wing2_altb_time 


A-9 


end 


end 

y,  Input  Matrix 

bl=[0  [(y_nominal/wingl_Tsi_time)*(pi/180)]  0; 
[l/wing2_vel_time3  0  0; 

0  [-(x_nominal/wingl_Tsi_time)* (pi/180)]  0; 

0  [(l/wingl_Tsi_time)]  0; 

0  0  0; 

0  0  Cl/(wing2_alta_time*wing2_altb_time)]] ; 

y,  Disturbance  Matrix 
b2=[l  0  0; 

0  0  0; 

0  nominal_vel* (pi/180)  0; 

0  0  0; 

0  0  0; 

0  0  - [l/(wing2_alta_time*wing2_altb_time)]] ; 

B=[bl  b2]  : 

C*[l  0  0  0  0  0; 

010000; 

001000; 

000100; 

0  0  0  0  1  0]; 

D=[0  0  0  0  0  0; 

000000; 
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000000; 

000000; 

0  0  0  0  0  0]; 


I=C0  0  0  0  0  0]»; 


*/,  Closed  Loop  plant 
G=pi/180; 

al=[0  -1  y_nominal*((q*S*delta_C_Y_W_y/(m*nominal_vel))  -  G*KYP*Ky/tpsi)  -(G/tpsi)*y_nomi 

(q*S/(m>i<nominal_vel))*(y_nominal*delta_C_Y_W_z)  0  0  y_nominal*KYI*G/tpsi  0]; 

a2=[-KXP*Kx/tv  -(l/tv)*(l  +  KXP*Kv)  (q*S/m)*delta_C_D_W_y  000  KXI/tv  0  0]; 

a3=[0  0  x_nominal*(G*KYP*Ky/tpsi  -  (q*S/(m*nominal_vel))*delta_C_Y_W_y)  ... 

G*(x_nominal/tpsi  -nominal.vel  +  x_nominal*KYP*Kpsi/tpsi)  ... 

-x_noniinal*(q*S/(m*nominal_vel))*delta_C_Y_W_z  0  0  -x_nominal*KYI*G/tpsi  0]; 

a4=[0  0  (-KYP*Ky/tpsi  +  (q*S/(m*nominal_vel*G))*delta_C_Y_W_y)  -(l/tpsi)*(l+KYP*Kpsi)  ... 

(q*S/(m*nominal_vel))>i'delta_C_Y_W_z/G  0  0  KYI/tpsi  0]; 

a5=[0  0000100  0]; 

a6=[0  0  (q*S/m)*delta_C_L_W_y  0  -(l/(tha*thb))*(KZP*Kz  +  1)  -(1/tha  +  1/thb)  0  0  KZI/(tha 
a7=C-Kx  -Kv  0  0  0  0  0  0  0] ; 
a8=[0  0  -Ky  -Kpsi  00000]; 
a9=[0  0  0  0  -Kz  0  0  0  0] ; 

Acl=  [al ;  a2 ;  a3 ;  a4 ;  a5 ;  a6 ;  a7 ;  a8 ;  a9]  ; 
eig(Acl) 

open_system  controller 
cd  ../../test/; 
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time=input( ^please  enter  time  for  simulation  in  increments  of  10  seconds  :0; 
if  time  ==  0; 
time  =  def.time; 
else 

time  =  time* 10; 

end 

y,  Set  initial  state  values 

y,xi  =  [0 ;  nominal.vel ;  nominal.vel ;  x.nominal ;  y.nominal ;  0 ;  0] 

y, ;  nominal  .heading ;  nominal.vel ;  nominal.heading] 

y,x2  =  [0 ; nominal. vel ;  0 ;  0 ;  0 ; nominal. alt ;  nominal.alt]  ; 

y.x3=  []; 

y.XI  =  Cxl;x2]; 

max_error=le-3; 

min_step=t ime/ 1000 ; 

max.step=time/1000; 

y*Ct  ,x,y]  =rk45(' controller’ , time, XI,  Cmax.error,min_step, max. step]  )  ; 
y#Ct,x,y]  =rk45(’ controller’  ,time) ; 
yiCt,x,y]=rk45(’cont .generic’  ,time)  ; 

[t,x,y]=adams(’ controller’ ,time) ; 

y»[t ,x,y]  =rk45(’ controller’ , time, XI,  Cmax.error,min.step,max.step] ) ; 
y,  velocities  of  lead  and  wing 

set(f igure, ’Position’ ,  [560*(hx-l) ,480-480*hy,560  460] ) ; 
subplot (6, 1,1) ,  plot(t,y(: ,3) , ’-’) ; 
subplot(6,l,l) ,  ylabeK’Vel  (ft)’); 
hold; 

subplot (6, 1,1) ,  plot(t,y(: ,5) , ’c — ’) ; 
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*/,  headings  of  lead  and  wing 
subplot (6, 1,2) ,  plot(t,y(: ,4) 
subplot(6,l,2) ,  ylabelC’Psi  (deg)’); 
hold; 

subplot (6, 1,2)  ,  plot(t,y( :  ,6)  ,  ’c— ’)  ; 

*/,  lead  altitude 

subplot (6 ,1,3),  plot (t ,y ( : , 12) , ’ - ’ ) ; 
hold 

'/.  wing  altitude 

subplot (6 ,1,3),  plot (t  ,y ( ; ,  13)  ,  ’ c—  ’ )  ; 
subplot  (6,1, 3) ,  ylabeK’alt  (ft)’); 

'/,  X  seperation 

subplot  (6,1,4),  plot  (t  ,y  ( : ,  1)  ,  ’  c“  ’ ) ; 
hold; 

subplot (6, 1,4) ,  line(t ,x_nominal.*ones(size(t))) ; 
subplot (6, 1,4),  ylabeK’x  (ft)’); 

*/,  y  seperation 

subplot (6, 1,5) ,  plot(t,y(:  ,2)  ,  ’— o’) ; 
hold; 

subplot (6, 1,5) ,  line(t,y_nominal.*ones(size(t))) ; 
subplot (6, 1,5),  ylabeK’y  (ft)’); 

'/,  2  seperation 

subplot (6, 1,6),  plot(t,y(;  ,11) ,’— c’)  ; 
hold 

subplot (6, 1,6) ,  plot(t,z_nominal.*ones(size(t))) ; 
subplot(6,l,6)  ,  ylabeK’z  (ft)’); 
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xlabeK’Time  (sec)’); 


hx=hx+l; 

last_y=size(y,l) ; 

set(figure, ’Position’ ,  [560*(hx-l) ,480-480*hy,560  460]); 
plotCyC: ,2) ,y(: ,1)) ; 
hold  on 

plot (y (1,2) ,y(l,l) , ’o’) ; 
hold  on 

plot(y(last_y,2) ,y(last_y,l) , ’x’) ; 
ylabeK’x  (ft)  ’) ; 
xlabeK’y  (ft)’); 

a=x; 

•/.  TRUE  CHANGE  IN  DRAG 
'/)  nominal 

G=(2*Cl*nominal_vel*b)/(pi*AR) ; 

W=G/(4*pi*bp)*(log((y_nominal“2  +  z_nominal~2  +  (b*mu)“2)/((y_nominal-bp)“2  + 
z_nominal~2  +  (b*mu)“2)) . . . 

-  log(((y_nominal+bp)''2  +  z_nominal“2  +  (b*mu) ~2)/(y_nominal“2  +... 
z_nominal''2  +  (b*mu)“2))); 

delta_C_D=Cl*W/nominal_vel ; 

D_w_nom=q*S*delta_C_D/m; 

'/,  final 

G=(2*Cl*coinmand_vel*b)/(pi*AR)  ; 

W=G/(4*pi*bp)*(log((y_coinmand"2  +  z_command“2  +  (b*mu)“2)/((y_command-bp)''2  + 
z_command'‘2  +  (b>t'mu)~2))  . . . 

-  log(((y_command+bp)~2  +  z_coinmand~2  +  (b*mu)“2)/(y_command~2  +... 
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z_conmiand“2  +  (b*mu)“2))); 

delta_C_D=Cl*W/coinmand_vel; 

D_w_com=q*S*delta_C_D/m; 

'/,  change  in  drag 
delta_diff=(D_w_com-D_w_nom) 


A, 3  Linear  Simulation  Model 

The  linear  model,  shown  in  Figure  A.l  is  the  State-space  representation  of  the 
system  described  in  section  4.1  and  described  by  Equations  (4.13)  -  (4.16).  The  com¬ 
mands  input  to  the  Lead  are:  commanded_velocity-nominaLvelocity,  commanded^heading- 
nominalJieading,  and  commanded_alt-nominal^lt. 

The  individual  blocks  that  make  up  the  linear  model  are  described.  The  Wing&kinematic&vortex 
block  is  the  state  space  system  described  by  Equations  (4.13)  -  (4.16).  The  command  pre¬ 
filter  block,  Figure  A. 2  contains  the  prefilter  discussed  in  Chapter  V.  The  FFC  controller 
is  made  up  of  the  linear  mixer  block  described  in  Figure  A. 5  and  the  PI  controller.  Figure 
A. 6  described  in  Chapters  II  and  IV.  The  Lead  aircraft.  Figure  A. 3  is  composed  of  the 
flight  control  system  detailed  in  Chapter  11.  Finally,  the  Perts+nominals  block.  Figure 
A.4  is  used  to  add  the  nominal  settings  back  to  the  perturbations  to  obtain  the  actual 
separation  distances. 


Figure  A.l  Linear  State  Space  Matlab  Simulation  Model 
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Figure  A. 3  Lead  Aircraft  Flight  Control  System  Model 
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Figure  A.4  Perturbation+Nominal  Block 
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Figure  A.5  Linear  Mixer  Simulation  Model 
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Figure  A.6  PI  Controller  Simulation  Model 
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A. 4  Non-Linear  Simulation  Model 

The  full  block  simulation  model  is  shown  in  Figures  A.7  -  A. 10.  The  prefilter,  linear 
mixer,  PI  controller,  and  lead  aircraft  blocks  are  the  same  as  those  used  in  the  linear 
model  above.  The  Kinematic  and  Horseshoe  Vortex  simulation  blocks  contain  both  the 
linear  models  and  the  non-lineax  models.  Figures  A.9  and  A.IO  respectively.  The  full 
nonlinear  model  was  discussed  in  section  4.5.  The  wing  aircraft,  Figure  A.8  is  the  same  as 
the  Lead  except  the  disturbances  are  added  to  the  respective  autopilots  they  affect. 


y3 

Figure  A.7  Matlab  Block  Simulation  Model 
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Figure  A. 8  Wing  Aircraft  Flight  Control  System  Model 


A-22 


Figure  A.9  Kinematic  Motion  Simulation  Model 
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Figure  A.IO  Horshoe  Vortex  Disturbance  Simulation  Model 
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In  this  paper  the  close  formation  flight  control  problem  is  addressed.  The  formation  consists  of  a  lead  and  wing  aircraft, 
where  the  wing  flies  in  close  formation  with  the  lead,  such  fliat  the  lead's  vortices  produce  aerodynamic  coupling  effects, 
and  a  reduction  in  the  formation's  drag  is  achieved.  A  controller,  i.e.,  a  formation-hold  autopilot  for  the  wing  aircraft,  is 
designed  such  that  the  formation's  geometry  is  maintained  in  the  face  of  lead  aircraft  maneuvers.  In  the  formation  flight 
control  system,  the  wing  and  lead  aircraft  dynamics  are  coupled  due  to  kinematic  effects,  and,  in  the  case  of  close 
formations,  additional  aerodynamic  coupling  effects  are  introduced.  In  the  paper  these  additional  aerodynamic  coupling 
effects  are  properly  modeled.  The  most  significant  aerodynamic  coupling  effect  introduced  by  close  formation  flight 
entails  the  coupling  of  the  lateral/directional  channel  into  the  altitude-hold  autopilot  channel.  It  is  shown  that  formation 
hold  autopilots  designed  ignorin  the  aerodynamic  coupling  effects  yield  satisfactory  performance  in  close  formation 
flight. 
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Close  Formation  Flight  Control,  Formation  Auto-pilot,  Formation  Autopilot,  F-16,  Flight 
Control,  Horseshoe  Vortex,  UAV 
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